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I.  INTRODUCTION 


This  report  describes  a  research  program  to  determine  the  sources 
of  scattering  in  Ti-dif fused  LiNbO^  waveguides,  to  eliminate  or  reduce  these 
sources  and  their  effects ,  and  to  evaluate  the  finished  waveguides  for  their 
potential  use  in  the  integrated  optical  spectrum  analyzer  (IOSA)  application/1^ 
The  motive  for  the  research  is  the  fact  that  light  scattered  in  the  plane  of 
the  waveguide  can  arrive  at  unintended  elements  in  the  IOSA  detector  array. 
Scattering  is  therefore  a  source  of  noise  that  limits  the  smallest  detectable 
information  signal,  or  equivalently,  the  dynamic  range  of  the  system.  Reduc¬ 
tion  in  waveguide  scattering,  especially  in-plane  waveguide  scattering,  is 
therefore  useful  in  improving  the  dynamic  range  of  systems  like  the  IOSA,  and 
in  helping  these  systems  to  attain  their  maximum  theoretical  potential. 

The  logical  plan  for  scattering  reduction  research  involves  the  use 
of  diagnostic  tools  for  "visualizing"  the  scattering  centers,  the  development 
of  waveguide  fabrication  or  processing  techniques  to  reduce  the  density  and/or 
the  scattering  cross  section  of  these  scattering  sources,  the  development  of 
waveguide  evaluation  procedures  for  verifying  the  quality  of  the  waveguides, 
and  a  theory  to  relate  the  measured  level  of  scattering  to  the  performance  of 
integrated  optical  systems,  such  as  the  IOSA,  that  will  employ  the  waveguides. 

We  report  only  limited  success  in  achieving  the  first  three  of  these 
tasks,  experimental  in  nature.  A  variety  of  standard  diagnostic  tools  have  been 
tested  for  their  ability  to  provide  information  about  the  size  and  density  of 
scattering  centers.  These  include  differential-interference-phase-contrast 
microscopy,  scanning  electron  microscopy,  e-beam  X-ray  analysis,  and  multiple- 
beam-interference  microscopy.  In  general,  we  find  direct  evidence  of  scatter¬ 
ing  centets  only  in  waveguides  that  have  been  purposely  fabricated  under  adverse 
conditions.  Conventional  fabrication  methods  produce  waveguides  that  have  no 
apparent  scattering  centers  as  determined  by  the  diagnostic  methods  listed  above. 

This  negative  result  has  motivated  us  to  look  for  scattering  centers 
by  indirect  methods.  The  measurement  of  in-plane  scattering  is  one  such  method. 
It  is  an  important  one  in  that  it  provides  the  ultimate  information  as  to  how 
the  waveguide  will  function  in  an  integrated  optical  signal  processing  applica¬ 
tion.  With  ln-plane  scattering  measurements,  we  can  at  least  predict  how  a 
waveguide  will  perform  in  a  specific  application,  even  if  we  can’t  identify 
the  actual  scattering  centers  at  work. 


In-plane  scattering  can  further  serve  as  a  diagnostic  tool  if  care¬ 
ful  models  of  waveguide  scattering  centers  are  constructed.  To  do  this,  we 
identify  those  steps  in  the  waveguide  fabrication  procedure  where  inhomogene¬ 
ities  are  likely  to  be  introduced  into  the  material.  We  estimate  or  otherwise 
determine  the  refractive  index  difference  associated  with  a  specific  inhomo¬ 
geneity,  and  then  we  determine  if  reasonable  values  of  scattering  center  size 
and  density  on  the  waveguide  surface  produce  a  calculated  level  of  in-plane 
scattering  consistent  with  measured  levels.  If  so,  the  scattering  mechanism 
being  considered  is  a  plausible  one  worthy  of  further  consideration.  This  may 
involve  modifying  the  waveguide  fabrication  process  to  effect  a  change  in  the 
candidate  scattering  mechanism,  and  then  checking  to  see  that  the  in-plane 
scattering  is  similarly  affected.  A  good  example  is  that  of  scattering  by  sur¬ 
face  roughness.  If  measured  levels  of  scattering  are  consistent  with  a  certain 
level  of  surface  roughness,  then  polishing  of  the  surface  should  result  in  a  de¬ 
creased  level  of  in-plane  scattering.  If  it  does  not,  volume  scattering  centers 
are  indicated,  and  a  different  set  of  fabrication  procedures  may  be  employed  to 
learn  the  source  of  these. 

This  indirect  approach  has  provided  the  main  thrust  of  our  work  in 
the  absence  of  useful  diagnostic  tools.  However,  even  this  approach  has  been 
frustrated  by  the  difficulties  encountered  in  taking  reliable  in-plane  scatter¬ 
ing  data.  The  measurements  invariably  contain  a  background  scattering  contri¬ 
bution  associated  with  input  and  output  coupling  from  the  waveguide.  This 
contribution  compares  with  and  in  some  cases  dominates  the  scattering  from 
the  waveguide  that  we  wish  to  measure.  As  a  result,  the  error  bars  on  the 
experiments  are  too  large  to  permit  us  to  draw  conclusions  regarding  the  effi¬ 
cacy  of  various  modifications  of  the  waveguide  fabrication  procedure.  As  far 
as  we  can  tell  from  our  experiments,  a  wide  range  of  Ti-film  thicknesses,  dif¬ 
fusion  times  and  temperatures  are  useful  in  producing  equally  good  waveguides. 

The  most  promising  interpretation  of  our  experimental  results  shows  that  the 
waveguides  are  better  than  required  to  meet  dynamic  range  specifications  for 
an  integrated  optical  spectrum  analyzer.  The  most  negative  interpretation  is 
that  5  dB  improvement  is  still  required. 

For  the  near  term  then,  the  question  as  to  whether  waveguide 
quality  is  sufficient  for  integrated-optical-systems  applications  will  have  to 
be  answered  by  the  fabrication  and  characterization  of  those  systems.  The 
present  research  will  be  useful  as  a  starting  point  to  those  who  may  continue 
it.  Additionally,  we  have  identified  two  promising  candidates  for  in-plane 
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scattering  and  eliminated  several  others.  The  two  promising  candidates  are 
surface  roughness  and  near-surface  inhomogeneities  associated  with  the  diffu¬ 
sion  process.  The  latter  may  be  the  lithium  titanate  compounds  identified  by 
Burns  et  al(2i?sing  Secondary  Ion  Mass  Spectroscopy  (SIMS).  The  relative  impor¬ 
tance  of  these  candidate  scattering  mechanisms  may  be  determined  by  performing 
SIMS,  Talystep  profilometry  and  in-plane  scattering  measurements  on  samples 
fabricated  using  a  variety  of  procedures.  This  was  in  fact  the  plan  of  the 
most  recent  research,  but  it  was  not  implemented  in  view  of  the  in-plane  scat¬ 
tering  experimental  problem.  Our  effort  was  reduced  to  in-plane  scattering 
measurements  on  existing  waveguides  to  determine  as  accurately  as  possible  the 
suitability  of  those  waveguides  for  an  integrated  optical  spectrum  analyzer. 
The  origin  of  scattering  centers  in  the  waveguides  was  not  addressed  experi¬ 
mentally. 

The  organization  of  this  report  is  as  follows.  In  Section  II,  we 

derive  the  IOSA  dynamic  range  dependence  on  generalized  scattering  parameters 
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such  as  E,  the  surface  density  in  cm  of  scattering  centers,  and  o(<}>)  the 
differential  scattering  cross  section.  Then  we  derive  expressions  for  o(<j>), 
in  terms  of  the  characteristic  size  and  refractive  index  of  the  scattering 
centers.  This  leads  us  to  the  consideration  of  three  distinct  scattering 
regimes:  Rayleigh  scattering  from  inhomogeneities  small  compared  to  a  wave¬ 
length,  Rayleigh-Gans-Debye  scattering  from  inhomogeneities  larger  than  a 
wavelength  but  with  a  small  refractive  index  difference  from  the  environment, 
and  Mie  scattering  from  inhomogeneities  larger  than  a  wavelength  and  having 
a  large  refractive  index  difference  from  the  environment. 

In  Section  III,  we  will  consider  a  variety  of  possible  scattering 
mechanisms  including  surface  roughness,  the  formation  of  lithium-titanate 
compounds,  the  formation  of  the  separated  phase  LiNbo0o,  and  nonuniform  dis- 

J  O 

tribution  of  diffused  Ti.  We  will  select  a  range  of  scattering  parameters 
associated  with  each  mechanism  and  use  these  to  calculate  theoretical  values 
of  IOSA  dynamic  range.  This  exercise  is  useful  in  that  it  tells  us  which 
scattering  mechanisms  are  likely  to  have  the  most  impact  on  IOSA  dynamic 
range,  and  which  mechanisms  are  likely  to  be  unimportant  even  if  they  are 
present.  Perhaps  the  most  significant  piece  of  information  generated  here 
is  the  likelihood  that  Rayleigh  scattering  is  not  a  factor  limiting  dynamic 
range  if  waveguide  attenuation  losses  are  1/4  dB/cm  or  less. 
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In  Section  IV  we  will  present  our  experimental  results  and  interpret 
them  in  relation  to  the  analysis  of  Sections  II  and  III.  We  will  begin  with 
results  showing  the  total  integrated  scattering  measured  as  a  function  of 
prism  separation.  These  results  suggest  the  importance  and  even  the  dominance 
of  prism-associated  scattering  to  the  total  scattering  generated  in  the  free 
waveguide.  We  will  describe  attempts  to  avoid  this  problem  using  end-fire 
coupling  techniques,  longer  waveguide  paths,  and  a  novel  optical  geometry  that 
distinguishes  between  waveguide  scattering  and  scattering  generated  at  the 
point  of  input  coupling.  Finally,  we  will  discuss  results  corrected  for  prism 
scattering  by  acquiring  data  at  two  different  coupling  prism  separations,  and 
then  subtracting  to  get  the  free-waveguide  contribution.  The  results  will  be 
associated  with  waveguides  formed  using  a  wide  variation  of  fabrication 
conditions  and  waveguides  fabricated  at  another  laboratory. 

In  Section  V  we  present  our  conclusions  as  to  likely  sources  of 
scattering,  best  methods  for  fabricating  waveguides  to  minimize  the  effects 
of  these  sources,  and  finally,  the  prospects  for  achieving  an  integrated- 
optical  spectrum  analyzer  having  the  desired  40  dB  dynamic  range. 
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II.  THE  EFFECTS  OF  OPTICAL  WAVEGUIDE  SCATTERING 
ON  IOSA  DYNAMIC  RANGE 

CHARACTERIZATION  OF  THE  SCATTERING  PHENOMENON 

The  effect  of  in-plane  waveguide  scattering  on  the  dynamic  range  of 
an  IOSA  depends  on  the  design  geometry  of  the  device  and  on  the  scattering 
parameters  of  the  waveguide.  The  scattering  effect  of  a  single  scattering 
center  may  be  characterized  by  a  differential  cross  section  a (<}>)  defined  as 
follows:  If  I  is  the  optical  power  per  unit  width  of  the  confined  beam  at 
the  location  of  the  scattering  center,  then 

AP  -  Io(4>)A<j>  (1) 

is  the  power  scattered  in  the  plane  of  the  waveguide  at  an  angle  cp  into  an 
angular  sector  A$.  The  total  scattered  power  is  obtained  by  summing  the  in¬ 
dividual  contributions  described  by  Eq.  (1).  This  is  most  readily  done  by 
taking  o(< J>)  to  be  a  statistically  averaged  cross  section  for  the  ensemble  of 

scattering  centers,  by  taking  the  scattering  centers  to  be  uniformly  distributed 
_2 

with  a  density  I [cm  ]  over  the  waveguide  surface,  and  by  integrating  Eq.  (1) 
over  the  area  spanned  by  the  waveguided  beam.  The  total  power  scattered  into 
the  sector  A<j>  at  an  angle  <J>  is  found  to  be 

Aptot  *  PELo(*)A<t>  ,  (2) 

where  L  is  the  path  length  and  P  is  the  power  in  the  unscattered  beam.  It  is 
assumed  that  P  is  not  significantly  depleted  over  the  path  length  by  either 
scattering  or  absorption  phenomena. 

A  measurmement  of  APtQt  versus  <f>  enables  us  to  determine  the  product 
£a(<J>).  With  this  information  and  knowledge  of  the  device  geometry,  we  can 
predict  the  dynamic  range  of  an  IOSA.  The  mathematical  analysis,  as  we  will 
present  it,  requires  the  specification  of  the  functional  form  of  a < ) •  This 
depends,  in  general,  on  the  size  and  refractive  index  of  the  scattering  centers. 
However,  certain  simplifications  result  in  the  special  cases  of  Rayleigh  scat¬ 
tering  (RS)  and  Rayleigh-Gans-Debye  scattering  (RGDS) .  For  the  latter,  we  will 
use 
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(3) 


o(t)  -  cr(0)  [1  +  (ka<fr)2]-1 

where  k  *  2ir/X,  X  is  the  wavelength  in  the  material,  a  is  the  characteristic 
scattering  center  size,  (assuming  a>>X  for  RGDS),  and  0  is  the  scattering  angle 

in  the  material,  given  in  radians.  The  Lorentzian  form  for  a($)  has  been  de- 

(3) 

rived  by  Marcuse,  and  we  will  later  show  its  appropriateness  using  a  some¬ 
what  different  analysis.  For  Rayleigh  scattering,  the  form  of  the  scattering 
cross  section  is  complicated  by  the  waveguide  geometry.  We  will  use  formulas 
derived  for  uniform  three-dimensional  space,  but  applied  only  in  the  plane  of 
the  waveguide.  They  are 

2 

cr(<j>)  *  o(0)  cos  <p,  (TE  polarization),  (4) 

cr(<j>)  ■  o(0)  ,  (TM  polarization). 

The  case  of  Mie  scattering  is  too  complicated  to  be  dealt  with  rigorously 
here.  Generally  the  angular  range  will  be  intermediate  to  the  small  angular 
range  for  RGD  scattering  and  the  large  range  for  Rayleigh  scattering.  This 
will  enable  us  to  qualitatively  discuss  the  case  of  Mie  scattering  in  relation 
to  our  findings  for  Rayleigh  and  RGD  scattering. 

CHARACTERIZATION  OF  THE  IOSA  GEOMETRY 


Spectrum  analyzer  dynamic  range  depends  on  design  parameters  as 

well  as  scattering  parameters.  Figure  1  shows  the  system  geometry  used  in 

(4) 

the  Westinghouse  IOSA.  The  actual  numerical  values  for  the  parameters 
listed  in  Fig.  1  are  shown  in  Table  I.  We  have  called  attention  to  four  dif¬ 
ferent  regions  of  the  spectrum  analyzer  where  separate  calculations  of  the 
scattered  power  incident  on  the  detector  array  are  to  be  made,  with  the  re¬ 
sults  added  to  obtain  the  total  scattered  power.  In  regions  C  and  D,  two 
beams  are  present,  called  1  and  0,  standing  respectively  for  the  diffracted 
(lst-order)  and  undiffracted  (0-order)  beams  obtained  after  passage  of  the 
guided  beam  through  the  surface  acoustic  wave  that  divides  regions  B  and  C 
of  the  IOSA.  The  scattering  contribution  from  each  of  these  beams  must  be 
considered  separately,  so  there  are  a  total  of  six  different  scattered  powers 
to  be  added  to  obtain  the  total  scattered  power  at  the  various  elements  of 
the  IOSA  detector  array.  We  refer  to  these  powers  as  P^,  where  the  subscript  i 
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TABLE  I.  PARAMETERS  AND  PARAMETER  VALUES  USED  IN 
IOSA  DYNAMIC  RANGE  ANALYSIS 


Parameter 


Description 


Numerical 

Values 


» 

I 


# 

n 

Sy 

Ay 

n 

k 

a 

o(0) 

£ 


length  of  waveguide 

collimating  lens  focal  length 

lens-SAW  separation 

Fourier-transform  lens  focal  length 

AO  diffraction  angle 

AO  diffraction  efficiency 

pixel  width 

pixel  spacing 

mode  index 

free-space  wavelength 
2nn/X0,  wave  vector  magnitude 
scattering  center  size 

peak  differential  scattering  cross  section 
scattering  center  surface  density 


72.4  mm 
27. 2  mm 
9  mm 

27.2  mm 
0.066  radians 

5% 

12  ym 
12  um 

2.2 

0.83  um 
16.65  um"l 
various 
various 
various 


is  AO,  BO,  CO,  Cl,  DO,  or  Dl,  depending  on  the  region  and  the  beam 
that  generates  the  scattering.  Some  simplification  of  the  calculation  is 
afforded  by  the  fact  that  subscripts  AO,  DO,  and  Dl  are  associated  with 
cylindrical  waves,  while  BO,  CO,  Cl  are  associated  with  plane  waves.  We 
consider  the  scattering  from  these  first. 

Scattered  Power  from  Plane-Wave  Portions 
of  the  Guided  Beam 

For  power  scattered  from  plane  waves  into  an  angular  sector  AiJ> 
located  at  $,  the  Fourier-transform  lens  in  Fig.  1  will  deliver  the  power  to 
a  location  removed  from  the  beam  axis  in  the  plane  of  the  detector  array  by  a 
distance  F<p .  The  power  will  be  distributed  over  a  width  FA<p .  We  will  define 
the  transverse  coordinate  in  the  array  plane  to  be  Y,  where  Y  *  0  is  located 
at  the  central  element  in  the  detector  array.  If  each  element  has  a  width 
6Y,  the  subtended  scattering  angle  is  ■»  6Y/F.  For  the  diffracted  beam  Cl, 
the  actual  scattering  angle  is  <t>  *  Y/F.  The  scattered  power  received  by  the 
detector  located  at  Y  is  therefore 

APC1(Y)  -  PC1  IV  o(Y/F) (6Y/F) ,  (5) 

where  L'  is  defined  in  Fig.  1  as  the  distance  of  the  acoustic  wave  center 
from  the  center  of  the  Fourier-transform  lens.  P  is  the  unscattered 
power  in  the  diffracted  beam  in  region  C.  If  we  set  pcl  “  p  and,  in 
general,  define  Pi  as  the  unscattered  power  for  all  of  the  beams  in  the  IOSA, 
where  i  “  AO,  BO,  CO,  Cl,  DO,  and  Dl,  the  relation  between  the  six 

values  of  P^  and  P  is  determined  from  the  acoustic-wave  diffraction  efficiency 
n  as  follows: 

P  AO  *  P  BO  "  P^n  > 

PC1  "  p  Dl  “  p  •  (6) 

PC0  -  p  DO  -  CP/o)  (1-n) . 


This  treatment  ignores  guided-wave  attenuation  associated  with  out-of-plane 
scattering  or  absorption. 

To  calculate  the  scattered  power  from  the  undiffracted  beam  in  region 
C,  we  use  the  third  of  these  relationships  along  with  the  fact  that  the  undif¬ 
fracted  guided  beam  propagates  at  an  angle  \l>  relative  to  the  optical  axis  defined 
by  the  direction  of  the  diffracted  beam.  As  a  result  light  which  impinges  at  a 
position  Y  in  the  detector  plane  is  scattered  through  an  angle  =  if)  -  Y/F.  The 
scattered  power  received  at  a  detector  centered  at  Y  and  associated  with  beam 
CO  is  therefore 


APC0(Y)  -  (P/n)  (l-n)ZL'aO|>-Y/F)  (6Y/F)  (7) 

The  scattered  power  from  the  collimated  beam  incident  on  the  acoustic  wave 
is  similarly  given  by 


APbo  *  (P/nWo^-Y/FMdY/F).  (8) 

The  total  scattered  power  from  the  collimated  beam  between  the  two  lenses  of 
the  IOSA  is  the  sum  of  Eqs.  (5),  (7),  and  (8).  Referring  to  this  scattered 
power  as  APB(,,  we  find 

APbc(Y)  -  PEL' (5Y/F)  [o(Y/F)  +  (2/n-l)o(*-Y/F) ] .  (9) 

Scattered  Power  From  Cylindrical-Wave  Portions 
Of  The  Guided  Beam 


The  analysis  is  more  complicated  for  beams  in  regions  A  and  D,  where 
the  wavefront  is  cylindrical  rather  than  planar.  (The  reader  will  have  noted 
that  we  are  considering  each  geodesic  lens  to  act  as  would  a  thin  lens  located 
at  the  geodesic  lens  center) .  Consider  first  the  scattered  power  contributed 
by  the  diffracted  beam  Dl.  Figure  2  shows  a  scattering  event  at  a  distance  z 
from  the  array  at  an  off-axis  position  y.  Integrating  the  contributions  from 
all  distances  z  and  positions  y,  we  find  the  scattered  power  distributed  over 
a  detector  element  centered  at  Y  to  be 
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Scattered 

Ray 


(10) 


F  «  VY*y*z) 

AP di  (Y)  -  S  /  dz  /  dyl(y,z)  /  d$  o(<J>)  , 

0  4>_(Y,y,z) 

where  the  upper  and  lower  limits  of  the  last  integrand  are  the  extremal  scat¬ 
tering  angles  subtended  by  the  detector  under  consideration.  For  most  values 
of  y  and  z  of  interest,  cr($)  will  not  vary  significantly  as  <p  varies  from  4>_ 
to  $+,  and  we  can  replace  the  last  integral  above  by  a($)  (<fi+-<j>  ),  where  <J> 
is  the  value  of  $  at  the  center  of  the  detector  element.  From  Fig.  2,  we  find 

4  «  tan-1 [ (Y-y) /z]  +  tan~1(y/z) 

=  tan_1[  (Y-y  +  5Y/2)/z]  (11) 

-  tan-1  [(Y-y  -  <5Y/2)/z] 

Even  this  simplification  is  not  enough  to  reduce  Eq.  (10)  to  a 
tractable  integral.  However  we  can  do  so  without  much  additional  loss  of 
accuracy  by  setting  y  »  0  in  Eq.  (11).  That  is,  we  take  the  scattering  angle 
<f>  and  the  subtended  angle  $+  -  <J>_  to  be  independent  of  the  off-axis  location 
of  the  scattering  event.  This  is  a  good  approximation  for  converging  beams 
subtending  a  low  numerical  aperture,  as  is  usually  the  case  in  planar  inte¬ 
grated  optical  applications.  We  find 

?  *  tan_1(Y/z)  (12) 

4>+-<f»_  *  tan-1  [  (Y  +  6Y/2)/z]  -  tan-1[(Y-<5Y/2)/z] .  (13) 

The  subtended  angle  can  be  further  simplified  to 

=  2tan-16Y/2z  ,  Y  -  0  (14) 

=  5Yz/(z2  +  Y2),  Y  t  0 

The  first  of  these  expressions  is  obtained  directly  by  setting  Y  -  0  in 

Eq.  (13).  The  second  is  obtained  by  rewriting  the  right  hand  side  of  Eq.  (13) 
-1  2  2  2 

as  tan  <5Yz/(z  +  Y  -6 y  /4),  and  then  replacing  the  inverse-tangent  function 


2 

by  its  argument.  Furthermore,  we  eliminate  the  6Y  / 4  term  in  comparison  to 

2  2  2  2  2  2 
Y  .  Since  Y  •  m  Sy  for  the  m  th  off-axis  detector,  where  m  >1,  Y  is  at 

least  four  times  larger  than  the  neglected  term. 

With  these  approximations,  the  scattered  power  from  the  diffracted 

beam  in  Region  D  is  given  by  the  integral 


F  j  2tan-1(<5Y/2z)l 

APn,  (Y)  -  I  /  dzf  dyl(y,z)a(tan-i  Y/z)'  or  /  ,  (15) 

0  -«  ! 6Yz/(z2  +  YZ) j 

where  the  top  (bottom)  term  in  curly  brackets  is  used  when  Y  equals  (does  not 
equal)  zero,  and  where  the  functional  form  of  a  is  that  of  either  Eq.  (3)  or 
(4),  depending  on  the  type  of  scattering  that  dominates.  The  integral  over 
y  is  readily  evaluated  to  give  the  power  in  the  optical  beam.  For  the  dif¬ 
fracted  beam  in  region  D, 


PD1  "  P  =  /  dy  I(y,z).  (16) 

— oo 

If  the  waveguide  is  lossless,  as  we  assume,  P  is  independent  of  z  because 
energy  is  conserved  at  each  distance  along  the  beam  axis.  For  lossy  wave¬ 
guides,  attenuation  mechanisms  effect  scattered  and  unscattered  light  similarly. 
Since  dynamic  range  is  calculated  as  the  ratio  of  scattered  to  unscattered 
power,  the  effects  of  losses  are  largely  self-cancelling.  This  shows,  in  a 
qualitative  way,  the  appropriateness  of  neglecting  waveguide  attenuation. 

After  inserting  Eq.  (16)  into  Eq.  (10),  and  incorporating  Eqs.  (3)  and  (4) 
for  the  various  forms  of  the  scattering  cross  section,  we  may  write  Eq.  (15) 
as 


APD1.(Y) 


=  PIc(O)  /  dz 
0 


[1 


+  k2a2(tan'1Y/z)2] 


2//  2  . 
z  f (z  + 


Y2) 


"W 


2tan~l(6Y/2z)'/ 
Wz/(z2  +  Y2)i 


(17) 


This  is  really  six  different  equations,  one  for  each  combination  of  terms  in 
the  two  sets  of  curly  brackets.  The  three  different  functions  in  the  first 
set  of  curly  brackets  are  associated  with  the  type  of  scattering  center  (RS 


or  RGDS)  and  the  polarization  of  guided  light  employed  (TE  or  TM,  but  TE  in 
the  IOSA  application) .  The  second  set  of  curly  brackets  contains  the  approxi¬ 
mate  angle  subtended  by  the  detector  element  centered  at  Y.  For  the  on-axis 
detector  at  Y  ■  0,  the  top  term  in  the  second  set  of  curly  brackets  applies, 
and  the  first  set  of  curly  brackets  contains  only  unity.  The  integral  over 
the  inverse  tangent  function  can  be  carried  out  exactly.  The  result  is 

APd1(0)  =  PZo(o)6Y[u  tan"1(l/u)  +  £n(l  +  u2)1/2],  (18) 

where  u  *  2F/6Y.  With  F  typically  several  centimeters  and  Y  typically  ten 
microns,  u  is  on  the  order  of  10^.  For  this  large  value,  Eq.  (18)  simplifies 
further  to 


ApDi(°)  =  PEo(0)SY  £n2F/5Y  (19) 

For  nonzero  values  of  Y,  Eq.  (17)  consists  of  three  different  in¬ 
tegrals  for  the  three  different  expressions  in  the  first  set  of  curly  brackets 

2  2 

multiplied  by  5yz/(z  +  Y  ).  The  two  integrals  for  Rayleigh  scattering  can 
be  calculated  analytically,  as  can  the  integral  for  RGD  scattering  if  we  re¬ 
place  tan  Y/z  by  Y/z  in  the  first  set  of  curly  brackets.  This  approxima¬ 
tion  is  invalid  at  small  z;  however,  the  integrand  vanishes  at  z  *  0  whether 
or  not  the  approximation  is  employed.  Consequently,  there  is  not  likely  to 
be  a  significant  error  in  the  integral  as  a  result. 

The  values  for  the  three  integrals  are  displayed  in  Table  II.  In 
two  of  the  three  expressions  there  is  a  second  term  that  is  generally  small 
in  comparison  to  the  first.  If  we  ignore  this  term  in  Table  II,  the  three  ex¬ 
pressions  acquire  the  same  form;  namely 

APD1(Y)  =  PIa(om  in[(F2  +  q2Y2)/q2Y2]1/2  ,  (20) 

where  q  *  ka>>l  in  the  case  of  RGD  scattering  and  q  ■  1  in  the  case  of  Rayleigh 
scattering,  independent  of  the  polarization  of  light  in  the  waveguide.  The 
expression  of  Eq.  (20)  diverges  as  Y-M).  However,  we  know  that  Eq.  (19)  should 
be  used  instead  of  Eq.  (20)  as  this  limit  is  approached.  It  is  useful  to 
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TABLE  II.  SCATTERED  POWER  FROM  THE  DIFFRACTED  BEAM  IN 
REGION  D,  FOR  DETECTOR  ELEMENTS  LOCATED  OFF- 
AXIS  AT  A  DISTANCE  Y 


Type  of  Differential  Scattering 

_ Scattering _ Cross  Section _ 

Rayleigh  Gans  a(0)  [1  +  k 2a2$2]  ^ 

Debye 


Rayleigh  2 

(TE  polarization)  o(0)  cosz  $ 


Rayleigh 

(TM  polarization)  a(0) 


APD1(Y)/PIa(0)6Y 

[k2a2/(k2a2-l)  ]£n[(F2+k2a2Y2)  /k2a2Y2]1/2 
-  [  1/  (k2a2-l)  ]  in  [(F2+Y2)  /Y2  ] 1/2 

£n[ (F2  +  Y2)/Y2)jL/2  -  (1/2)F2/(F2+Y2) 

£n[(F2  +  Y2)/Y2]1/2 


2  2 

phenomenologically  combine  the  two  results  by  replacing  q  Y  in  the  denominator 
of  Eq.  (20)  by  q^2  +  6y2/4.  The  added  term  removes  the  divergence  at  Y  *  0 
and  provides  the  correct  expression  as  that  limit  is  approached.  At  the  same 
time,  it  does  not  produce  a  significant  error  at  nonzero  values  of  Y  of  Interest, 
namely  Y  *  m6y,  where  m  *  ±1,  ±2,  ±3...±M  is  the  number  of  the  detector  element 
in  the  array  relative  to  the  on-axis  detector  element,  having  m  *  0.  In  summary 
the  scattered  power  at  detector  m  from  the  diffracted  beam  in  region  D  is  approxi¬ 
mately  given  by 

APD1(m)  -  PZa(0)<5YJln[(F2  +  q2m26Y2)/(5Y2/4  +  q2m25Y2) ]1/2.  (21) 

The  scattered  power  from  the  undiffracted  beam  in  region  D  is  obtained  from 
this  expression  by  replacing  m  by  Ftanty/6y  -  m  *  mo~m.  This  accounts  for  the 
propagation  of  the  undiffracted  beam  at  an  angle  relative  to  the  propagation 
direction  of  the  diffracted  beam.  We  also  replace  P  by  (P/n)(l-n)  to  account 
for  the  power  variation  in  the  diffracted  and  undiffracted  beams.  The  result 
is 

APD0(m)  -  (P/n)(l-n)Io(0)6Y  £n{  [F2  +  q2(mo-m)  2<5Y2] 

/ [6Y2/4  +  q2(mo-m)26Y2]}1/2  (22) 

The  undiffracted  beam  in  region  A  is  also  accounted  for  by  an  ex¬ 
pression  of  this  type.  This  is  not  Immediately  obvious  because  the  scattering 
in  this  region  is  produced  by  a  diverging  rather  than  a  converging  beam.  How¬ 
ever  one  effect  of  the  collimating  and  Fourier-transform  lenses  in  the  IOSA 
device  is  to  produce  in  the  detector  plane  an  image  of  the  light  pattern  in 
the  laser  plane.  The  effect  of  scattering  centers  in  region  I  is  to  produce 
cylindrical  wavefronts  that  may  be  regarded  as  emanating  from  virtual  sources 
in  the  laser  plane.  The  distribution  of  the  virtual  intensity  associated  with 
these  sources  is  calculated  by  propagating  the  scattered  wavefronts  backward 
through  the  scattering  centers  to  the  laser  plane.  The  virtual  image  that 
results  from  the  scattering  of  a  diverging  spherical  wave  is  therefore  identi¬ 
cal  to  the  real  image  that  results  from  the  scattering  of  a  converging  beam. 

The  IOSA  collimating  and  Fourier-transform  lenses,  in  producing  a  real  image 
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of  the  virtual  wavefront  in  the  laser  plane,  simply  replicate  the  scattered- 
intensity  distribution  from  the  undiffracted  beam  in  Region  D.  The  only 
difference  is  that  associated  with  the  variation  of  optical  power  in  the  two 
regions .  We  have 


APA0^  *  PD0(m)/(l-n) .  v 

The  Dynamic-Range  Formula 

The  total  scattered  power  at  detector  m  is  obtained  by  adding  the 
results  of  Eqs.  (9),  (21),  (22),  and  (23).  Neglecting  polarization  correc¬ 
tions  Eq.  (9)  can  be  written  in  the  form 


APBC(m)  *  PZL’(6Y/F)ct(0)  1 

{[1  +  q,2m26Y2/F2]_1  +  (2/n-l)  [1  +  q ' 2 (in^-m) 26Y2/F2 ]-1} 

where  q'  *  0  for  Rayleigh  scattering,  and  q'  *  k<X  for  RGD  scattering.  The 

2 

geometry  of  the  IOSA  is  such  that  the  coefficients  of  q'  are  much  smaller 
than  unity  for  all  detectors  m.  Consequently,  not  much  error  is  made  by 
replacing  q'  by  q  in  Eq.  (24),  where  q  *  ka  for  RGD  scattering  and  q  *  1 
for  Rayleigh  scattering.  This  enables  us  to  write  the  total  scattered 
power  in  a  unified  form,  valid  for  both  Rayleigh  and  RGD  scattering;  namely, 


APtot<m)/Pma(°)  - 

(L7F){[l+q2m2SY2/F2]-1  +  (2/n-l)  [l+q2(mQ-m) 2<5Y2/F2]-1} 
+in[(l+q2m26Y2/F2)/(q2m25Y2/F2  +  6Y2/4F2)]1/2 
+(2/n-l)An{(l+q2(mo-m)25Y2/F2]/[q2(-irj-m)26Y2/F2  +  SY2/4F2]  }1/2 


(25) 
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The  top  term  on  the  right  hand  side  is  the  contribution  from  the  plane  waves 
in  regions  B  and  C,  the  second  term  is  the  contribution  from  the  focused 
diffracted  wave  in  region  D,  and  the  third  term  is  the  contribution  from  the 
undiffracted  cylindrical  waves  in  regions  A  and  D. 

The  only  parameter  in  Eq.  (25)  that  relates  to  the  type  of  scatter¬ 
ing  mechanism  is  q,  where  q  =  1  for  Rayleigh  scattering,  and  q  *  ka  for  RGD 
scattering.  For  these  two  cases  at  least,  q  varies  approximately  as  the  in¬ 
verse  of  the  scattering  angle.  We  will  assume  that  this  situation  holds  for 
the  general  case  of  Mie  scattering.  That  will  enable  us  to  apply  Eq.  (24) 
to  any  scattering  mechanism  by  suitably  choosing  an  appropriate  value  of  q. 

Table  III  shows  the  design  parameters  for  an  IOSA  similar  to  the 

(4 ) 

one  reported  by  Mergerian  et  al.  Using  these  parameters  we  can  calculate 
the  relative  values  of  the  three  contributions  to  the  scattered  power  in 
Eq.  (25)  as  a  function  of  q,  taking  m  =  4.  This  choice  of  m  is  made  because 
IOSA  specifications  require  the  achievement  of  40  dB  dynamic  range  or  better 
at  the  4th  off-axis  pixel. 

TABLE  III.  DESIGN  PARAMETERS  FOR  AN  IOSA 


F  27. 2  mm 

L  18.0  mm 

6Y  12.0  um 

0.066  radians 
n  5% 

X  0.83  pm 


Lens  Focal  Lengths 
Lens-Lens  Separation 
Detector  Aperture  and  Spacing 
Diffraction  Angle 
Diffraction  Efficiency 
Free  Space  Optical  Wavelength 
Refractive  Index 


The  results  of  the  calculation  are  shown  in  Fig.  3.  At  small  values 
of  q,  the  dominant  contribution  comes  from  the  undiffracted  beam  in  regions 

A  and  D,  labeled  I  in  Fig.  3.  The  power  in  this  beam  is  about  20  times 
larger  than  in  the  diffracted  beam,  resulting  in  greater  scattered  power. 
Additionally  small  values  of  q  are  associated  with  large  scattering  angles, 
so  that  a  significant  portion  of  scattered  light  can  reach  the  detector  array. 

For  large  values  of  q,  the  scattering  angles  are  small,  and  light 
from  the  undiffracted  beam  does  not  reach  the  detector  array.  In  this  case 
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Relative  scattered  power  at  the  fourth  pixel  vs  q. 

I:  contribution  from  the  undiffracted  beam  in  regions 
A  and  D.  II:  contribution  from  the  diffracted  beam 
in  region  D.  Ill:  contribution  from  the  plane  waves 
in  regions  B  and  C. 


the  dominant  contribution  to  scattering-generated  noise  is  the  focuaed- 
diffracted  beam  in  region  0,  labeled  II  in  Fig.  3.  The  plane  wave  contribu¬ 
tion  from  regions  B  and  C  of  the  IOSA  is  labeled  III  in  Fig.  3  and  is  inter¬ 
mediate  to  the  other  contributions  at  most  values  of  q. 

IOSA  dynamic  range  is  computed  from  the  ratio  of  scattered  power 
captured  by  the  4th  detector  element  to  the  total  power  captured  by  the 
on-axis  detector.  Because  the  detector  elements  in  the  IOSA  array  are  large 
compared  to  the  focal  spot  size  of  the  beam,  the  power  illuminating  the  on- 
axis  detector  is  essentially  P.  The  corresponding  dynamic  range  is 
10  log[APtot(A)/P] .  This  may  be  calculated  from  Eq.  (25)  when  the  values  of 
Zcr(O)  and  q  are  specified.  The  use  of  Eq.  (25),  although  it  is  an  approxi¬ 
mation,  is  expected  to  provide  a  level  of  accuracy  commensurate  with  the 
accuracy  of  the  measurements  used  to  obtain  the  scattering  parameters. 
Furthermore  the  equation  is  simple  enough  to  program  on  most  common  program¬ 
mable  hand  calculators. 


CALCULATIONS  OF  SCATTERING  CROSS  SECTIONS  AND  DYNAMIC-RANGE 
PERFORMANCE  FOR  VARIOUS  TYPES  OF  SCATTERING  CENTERS 


While  values  of  Eo(0)  and  q  are  measured  experimentally,  it  is 
useful  to  calculate  the  dependence  of  o(0)  on  the  size  and  refractive-index 
inhomogeneity  of  the  scattering  center.  This  facilitates  the  use  of  the 
data  to  draw  conclusions  about  the  source  of  the  scattering.  We  consider 
three  different  situations:  (1)  Rayleigh  Scattering  (RS)  from  particles  of 
characteristic  size  &  small  compared  to  the  wavelength  (a<<\) ,  and  real  index 
of  refraction  different  from  that  of  the  immediate  environment  by  5n; 

(2)  Rayleigh  Gans  Debye  Scattering  (RGDS)  from  particles  large  compared  to 
the  wavelength  (CL»\)  and  refractive  index  such  that  Snka<<l;  (3)  Mie 
Scattering  (MS),  for  all  other  situations.  The  analyses  employed  to  calcu¬ 
late  o(0)  in  terms  of  Sn  and  a.  are  not  rigorous  but  are  supported  by  attrac¬ 
tive  physical  arguments.  Our  objective  is  to  obtain  the  correct  dependences 
on  the  parameters  involved  without  concern  for  exact  numerical  factors  or  for 
unnecessary  complications  associated  with  the  waveguide  geometry. 
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Figure  4  shows  a  schematic  drawing  of  Rayleigh  (dipole)  scattering 
centers  distributed  throughout  a  waveguide  layer  of  depth  D.  If  we  ignore 
the  effects  introduced  by  the  waveguide  geometry,  the  differential  scatter¬ 
ing  cross  section  is^ 

o' (a)  -  [k4a6(m2-l)2/(m2+2)2]cos2a,  (26) 

where  n/2-a  is  the  angle  between  the  scattering  direction  and  the  incident 
polarization,  and  m  ■  1  +  6n/n.  The  terms  containing  m  may  be  replaced  by 
the  approximation  (4/9) 6n  / n  .  While  this  approximation  is  obtained  in  the 
limit  of  small  <5n,  it  is  fairly  accurate  even  when  the  scattering  centers 
are  voids  having  n  +  <5n  =  1.0  and,  for  the  example  of  LiNbO^,  n*2.2.  The 
angle  a  in  Eq.  (26)  is  the  same  as  the  scattering  angle  4>  for  TE  polarized 
light  scattered  in  the  plane  of  the  waveguide,  while  a  *  0  for  TM  polarized 
light  scattered  in  plane.  This  has  been  used  in  writing  Eq.  (4);  however, 
o(0)  of  Eq.  (4)  can  not  be  identified  with  the  coefficient  o'(0)  of  cos  a 
in  Eq.  (26).  The  latter  coefficient  is  associated  with  scattering  into  a 
solid  angle  AO  in  three  dimensional  space,  while  o(0)  is  associated  with 
scattering  into  an  angular  sector  A$  in  the  plane  of  the  waveguide. 

The  relation  between  the  two  cross  sections  can  be  obtained  by 
deriving  the  power  scattered  into  a  waveguide  detector  by  two  equivalent 
methods,  one  utilizing  the  in-plane  cross  section  o( d>)  and  the  other  utiliz¬ 
ing  the  three-dimensional  cross  section  o' (a).  For  the  former,  the  inter¬ 
cepted  power  is  given  by  Eq.  (1).  For  the  latter,  the  intercepted  power 
generated  from  a  single  scattering  center  is 

AP  =  Zo'(a)Aft  ,  (27) 

where  I  is  the  power-per-unit-area  illuminating  the  scattering  center  and 
Aft  is  the  effective  solid  angle  subtended  by  the  detector.  For  a  guided- 
wave  power-per-unit  width  I  and  a  confinement  depth  D,  I  =  I/D.  In  calcu¬ 
lating  the  effective  solid  angle  subtended  by  a  detector,  Aft,  we  consider 
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only  scattered  rays  which  are  bounded  in-plane  by  the  detector  aperture  A# 
and  which  are  simultaneously  trapped  in  the  waveguide.  Trapped  rays  are 
generated  from  those  scattered  within  the  acceptance  angle  of  the  waveguide, 
A0.  This  angle  is  identified  in  Fig.  4,  and  is  found,  using  Snell's  Law, 
to  be 

A0  =  2(2An/n)1/2  ,  (28) 

where  An  is  the  index  change  associated  with  the  waveguide  layer  and  n>>An 
is  the  substrate  index.  The  effective  solid  angle  subtended  by  the  wave¬ 
guide  detector  elements  is  fA$A0.  The  factor  f  is  added  to  describe  the 
fact  that  only  a  fraction  of  rays  scattered  within  the  acceptance  angle  of 
the  waveguide  are  trapped  within  waveguide  modes. 

Using  these  results  in  Eq.  (27),  we  find  the  detected  power 

to  be 

AP  -  (I/D)  o'  (ot)  2f(2An/n)1/2A$  .  (29) 

For  power  scattered  in  the  plane  of  the  waveguide,  o' (a)  is  replaced  by 
2 

either  o'(0)cos  $  or  o'(0),  depending  on  whether  the  guided-wave  modes  are 
TE  or  TM.  By  comparison  of  this  result  with  Eqs.  (1)  and  (4),  we  make  the 
identification 


c( 0)  =  o' (0)  (2f/D)(2An/n)1/2  (30) 

This  expression  has  the  correct  dimension  of  [length]  \ 

The  dynamic  range  depends  on  the  product  Eo(0),  where  E  is  the 
number  of  scattering  centers  per  unit  area.  if  p  is  the  volume  density  of 
scatterers  in  the  waveguide,  then  E  ■  pD  and 

Eo(0)  -  po'(0)  (2f )  (2An/n)1/2  .  (31) 

It  is  further  useful  to  note  that 

o'  (0)  -  (3/8ir)otot  ,  (32) 
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where  0  is  the  total  scattering  cross  section  obtained  by  integrating 
0,(°)  over  all  solid  angles.  1  Substituting  this  result  into  Eq.  (31), 
we  obtain 

Zo(0)  -  (potot)(3/8TT)(2f)(2An/n)1/2  •  (33) 


Since  only  a  small  fraction  of  all  scattered  rays  are  retrapped  in  guided 
modes,  P°tot  ®ay  be  identified  with  the  waveguide  attenuation  coefficient 
associated  with  out-of-plane  scattering.  This  is  the  total  attenuation  co¬ 
efficient  in  waveguides  for  which  absorption  is  negligible,  and  for  which 
Rayleigh  scattering  centers  dominate.  Assuming  this  to  be  the  case  for  Ti- 
dif fused  LiNbO^  waveguides  used  in  IOSA  development,  Eq.  (33)  gives  us  the 
capability  of  calculating  IOSA  dynamic  range  performance  once  the  attenua¬ 
tion  coefficient  is  known.  Figure  5  shows  the  results  of  such  a  calcula¬ 
tion  for  the  case  f  *  1/2  and  An  =  0.01.  The  IOSA  dynamic  range  specifica¬ 
tion  of  40  dB  at  the  fourth  pixel  is  achieved  if  the  attenuation  coefficient 
is  less  than  0.25  dB/cm.  This  level  of  attenuation  should  be  readily  attained 
in  Ti:LiNb0.j  waveguides. ^  Consequently  Rayleigh  scattering  from  small 
inhomogeneities  should  not  be  an  important  factor  limiting  IOSA  dynamic- 
range  performance  in  systems  that  employ  TiiLiNbO^  waveguides. 

The  Case  of  Rayleigh-Gans-Debye  Scattering 

We  now  consider  scattering  centers  that  are  much  larger  than  a  wave¬ 
length,  but  which  are  optically  soft.  By  this  we  mean  that  the  field  inside 
the  scattering  center  is  essentially  the  same  as  the  incident  field,  and  the 
effect  of  the  scattering  center  is  simply  to  perturb  the  wavefront  of  the 
incident  light  by  an  amount  small  compared  to  a  wavelength.  These  require¬ 
ments  are  satisfied  when  q  ■  ka>>l,  and  ka5n<<l,  where  6n  is  the  index  per¬ 
turbation  associated  with  the  scattering  center. 

Figure  6  defines  the  geometry  of  the  problem  under  consideration. 

If  AQexp(ikz)  is  the  unperturbed  wavefront,  the  wavefront  after  passing 
through  the  scattering  center  is  AQexp(ikz)  +  Apgrt,  where  Apert  is  the 
perturbed  field.  If  the  scattering  center  has  a  characteristic  size  a,  then 
the  perturbed  field  just  beyond  the  scattering  center  is 


Range,  dB 


Apert  “  ikoa5ng  A0exp[ik(zo  t  a/2)],  in  yo-a/2<y<yQ+  a/2 


*  0,  elsewhere 


(34) 


In  our  notation  k  •  k  n  ,  where  k  ■  2ir/A  ,  A  is  the  free  space  wavelength, 

o  8  o  o  o 

n  Is  the  node  Index  of  the  guided  wave,  and  4n  is  the  change  in  mode  index 
8  8 
associated  with  the  material  index  perturbation  6n.  Equation  (34)  is  derived 

for  a  square-shaped  inhomogeneity  to  simplify  the  analysis  without  sacrificing 

its  important  features.  The  coordinate  pair  (yQ,zQ)  describes  the  center  of 

the  scattering  site. 

In  the  near  zone,  the  perturbed  field  of  Eq.  (34)  propagates  as  a 
plane  wave, but  in  the  far  zone,  diffraction  causes  the  wavefronts  to  become 
circles  centered  at  (yo,zQ).  The  form  of  the  perturbed  field  in  the  far  zone 
is  found  from  a  straightforward  application  of  Fraunhofer  diffraction  theory 
tobe<7> 


A(R,4>)  -  AoeikR(iAz)1/2ikoa25ngsinc(a<j,/A) 


e-ikyosin<t>  +  ikzo(l-cos<f) , 

where  R  is  the  radial  distance  from  the  scatterer,  z  is  the  distance  along 
the  axis  of  propagation,  and  $  is  the  angle  of  observation  relative  to  the 
direction  of  the  unperturbed  beam. 

For  an  ensemble  of  N  scattering  centers  the  total  scattered  field 
is  the  sum  of  N  expressions  like  that  of  Eq.  (35).  These  expressions  are 
identical  except  for  the  replacement  of  (yQ,zQ)  by  (y^.z^),  the  center  co¬ 
ordinates  of  the  ith  scattering  site.  The  result  of  squaring  the  magnitude 
of  the  total  scattered  field  is  a  randomly  fluctuating  intensity  distribution 
commonly  referred  to  as  a  speckle  pattern.  Because  the  longitudinal  dis¬ 
tances  between  the  scattering  centers  can  exceed  many  wavelengths,  the  speckle 

pattern  is  a  high-contrast  pattern  with  intensity  minima  near  zero  and  inten- 

(8) 

sity  maxima  which  can  be  several  times  larger  than  the  average  intensity.  The 
average  intensity  is  the  same  as  would  be  calculated  by  summing  the  intensi¬ 
ties  of  the  light  scattered  from  individual  scattering  sites,  and  is 

I  *  I  N(k  a2fin  )2(Xz)-1sinc2(a$/\) .  (36) 

scat  oo  8 
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The  angular  variation  of  the  scattered  intensity  with  <j>  is  contained  in  the 
2 

sine  (a$/X)  term.  However  this  function  itself  has  structure  that  results 

from  our  unphysical  assumption  of  square  scattering  centers  of  uniform  size. 

(3) 

Following  a  more  general  statistical  analysis  by  Marcuse,  the  formula 


I  ,  -  41  N(k  a26n  )(Xz)_1  L(ka*), 
scat  o  o  g 


L(ka$)  -  [1  +  (ka*)2]-1 


(37) 


is  derived.  Using  this  formula,  the  scattered  power  received  at  a  detector 
of  width  6Y  from  a  remote  scattering  area  of  length  L  and  width  W  is 


AP  •  I  fc6Y,  or 
scat  ’ 


AP  -  PSL4(koa26ng)2  X-1  A<f>L(ka<|>)  (38) 

We  have  used  N  ■  £LW,  P  *  IqW  and  A$  »  (<SY/z)  to  facilitate  comparison  of 
this  result  with  Eq.  (2).  The  comparison  enables  us  to  identify  the  differ¬ 
ential  scattering  cross  section  for  RGD  scattering  as 


o($)  ■  (2/ it)  n  k3a^5n2L(ka4) .  (39) 

go  g 

This  is  of  the  form  given  in  Eq.  (3),  with  o(0)  *  (2/ir)n  k3a^6n2. 

B  ®  B 

IOSA  dynamic  range  depends  on  the  product  Zo(0).  The  worst  case 

2 

situation  is  that  of  contiguous  scattering  centers,  for  which  I  *  1/a  ,  and 


Eo(0)  -  (2/ir)n  k3a26n2  .  (40) 

go  g 

Figure  7  shows  the  functional  relationship  between  d  and  An  that  will  re¬ 
sult  in  an  IOSA  dynamic  range  of  40  dB.  For  scattering  centers  between 
1  and  100  pm  in  size,  5ng  values  in  the  range  from  10  ^  to  2  x  10  or 
smaller,  are  indicated.  The  scattering  index  inhomogeneity  An,  in  other 
words,  should  perturb  the  guided  wave  mode  index  by  an  amount  three-to-four 
orders  of  magnitude  less  than  the  waveguide  index  change  An.  Just  what  values 
of  An  can  be  tolerated  depends  on  the  depth  of  the  scattering  center  beneath 
the  waveguide  surface,  its  thickness  measured  normal  to  the  surface,  and  the 
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modal  structure  of  the  guided  wave.  Thus  a  separate  analysis  is  required 
for  each  different  scattering  mechanism.  The  dependence  on  the  waveguide 
modal  parameters  suggests  the  possibility  of  reducing  in-plane  scattering 
by  appropriately  modifying  waveguide  fabrication  procedures  which  influence 
these  parameters. 

Calculation  of  fing  in  Terms  of  Scattering  Parameters 

The  perturbed  mode  index  depends  on  the  actual  index  inhomogeneity 
associated  with  the  scattering  center,  fin,  and  on  the  extent  to  which  the 
scattering  center  experiences  the  modal  field.  An  RGD  scattering  center  may 
be  characterized  as  having  a  characteristic  size  a  measured  in  the  plane  of 
the  waveguide  and  a  characteristic  thickness  zq  measured  normal  to  the  wave¬ 
guide.  If  zq<<D,  the  thickness  of  the  waveguide  layer,  the  scattering  center 
may  be  described  as  a  surface  scattering  center.  If  zq  “  a>>D,  it  may  be 
described  as  a  volume  scattering  center.  The  former  samples  only  the  surface 
modal  field,  and  may  be  reduced  in  Importance  by  generating  waveguides  having 
a  low  surface  field.  The  latter  samples  the  entire  mode,  and  may  only  be  re¬ 
duced  in  importance  by  reducing  fin  or  the  density  of  scattering  centers,  Z. 

These  considerations  are  made  quantitative  by  the  approximate 
relationship 


<Eq(z) | 5n(z) | Eo(z)>/<Eo(z) | Eq(z)> 


(41) 


where  Eq(z)  is  the  amplitude  profile  of  the  TEq  electric  field,  assumed  to 
be  employed  in  an  IOSA.  The  refractive  index  is 


I  5n(z)  *  fin 

-  0 


0<z<z 


elsewhere 


(42) 


where  z  *  0  is  the  waveguide  surface.  Utilizing  the  forms  for  Eq(z)  derived 
by  Marcuse  for  a  strongly  asymmetric  step  index  waveguide  having  a  small 
surface  index  change  An  and  an  air  superstrate,  we  find  for  zq<<D, 


6 tig  -  [4n/(n2-l)J  (k^fin)  S(Ang) 

S(An  )  -  (An-An  )(2n  An  ) 1/2  [1  +  k  D(2n  An  )1/2]“1  .  (43) 

8  8  8  ®  8 


30 


In  this  result,  S  is  a  function  proportional  to  the  modal  intensity 

at  the  surface.  It  uses  the  parameter  An  -  n  -  n,  where  n  is  the  substrate 

&  8 

index.  For  the  IOSA  application,  suitable  waveguides  are  those  for  which  An^ 

is  small  enough  to  preserve  single-mode  characteristics,  yet  large  enough  to 

(9) 

reduce  propagation  losses  associated  with  geodesic  lens  curvature.  We 

take  as  representative  n  ■  0.01,  An  ■  0.003,  and  D  *  2.15.  For  a  free  space 

®  -4 

wavelength  XQ  «  0.83  pm  and  n  =  2.2,  we  find  SCAn^)  =  2.80  x  10  .  The  mode 

index  perturbation  is  then 


<$n 


g 


6.42  x  10  ^  <5n  z 


(44) 


where  z  is  given  in  microns.  To  achieve  5n 
o  g 


,-6 


-3 

3.11  x  10  pm. 


2  x  10  ,  a  figure  for  which 

IOSA  dynamic  range  is  likely  to  be  above  40  dB,  we  require  <5n  zq 
In  Section  III  we  will  estimate  the  importance  of  some  real  scattering  mechanisms 
in  view  of  this  requirement. 

Note  that  for  zq>D,  Eq.  (43)  has  the  simple  result 


<5n  «  <5n  (45) 
g 

In  this  case,  the  actual  scattering  center  index  inhomogeneity  must  be  made 
less  than  2  x  10  to  achieve  a  high  probability  of  attaining  dynamic  range 
values  in  excess  of  40  dB. 


The  Case  of  Mie  Scattering 

For  present  purposes,  Mie  scattering  is  taken  to  describe  scatter¬ 
ing  from  index  inhomogeneities  much  larger  than  a  wavelength,  and  having  values 
of  6n  large  enough  to  justify  the  use  of  geometrical-optics  principles  in 
analyzing  their  scattering  effects.  If  the  typical  radius  of  curvature  of  a 
Mie  scattering  center  in  the  vicinity  of  the  guided  beam  is  a,  an  approximate 
application  of  the  lens-makers  formula  suggests  that  the  refracted  beam  will 
focus  at  a  distance  R  a  a/26n.  At  distances  greater  than  R,  the  refracted 
beam  diverges  with  a  half  angle  9  *  tan_1(a/R)  »  tan-1  2<$n. 

Figure  8  shows  a  highly  schematic  representation  of  the  event,  in 
which  the  diverging  beam  is  represented  as  a  cone  of  light.  The  top  portion 
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Waveguide 


of  the  light  cone  is  reflected  from  the  surface  of  the  waveguide  into  the  sub¬ 
strate,  doubling  the  average  intensity  in  that  region.  At  a  distance  r>>R,a 

from  the  scattering  center,  the  area  subtended  by  the  light  cone  is  a  semi- 

2  2 

circle  of  radius  rtan0  =  2r6n.  The  area  of  the  semicircle  is  A  *  2irr  <Sn  . 

If  the  power  scattered  by  the  imperfection  is  equal  to  the  power  Intercepted 
by  it,  then  AP  »  Id,  and  the  average  power  per  unit  area  in  the  scattered 
beam  is 

APtot/A  «  Id/2irr26n2  .  (46) 

The  effective  area  spanned  by  a  waveguide  detector  at  a  distance  r  is 
2 

A'  *  r  Aft,  where  Aft  is  the  effective  solid  angle  given  previously  as 

1/2 

fA$A0,  where  A9  -  2(2An/n)  is  the  acceptance  angle  of  the  waveguide,  and 

f  is  the  fraction  of  light  scattered  within  this  angle  that  is  trapped  in 

waveguide  modes  (dashed  rays  in  Fig.  8).  As  long  as  the  detector  is  within  the 

light  cone,  the  power  intercepted  by  it  is  AP  »  (AP  /A)A’ ,  or 

tot 

AP  »  (l/2ir)  2f(2An/n)1/2  (Id/ 6n2) A<f>.  (47) 

Outside  the  light  cone,  the  intercepted  power  is  zero.  By  comparison  with 
Eq.  (2)  the  differential  scattering  cross  section  is  therefore 

o($)  *  AP/IA<j)  *  (l/2ir)  (2An/n)^//2  (df/6n2),  in-  -tan  ^26n<<)><tan  ^26n 

*  0»  elsewhere.  (48) 

Calculations  of  IOSA  dynamic  range  require  knowledge  of  Eo(o)  and 
q,  where,  as  indicated  earlier,  q  is  correlated  with  the  inverse  of  the  scat¬ 
tering  half  angle.  For  the  example  of  Mie  scattering,  we  will  take  q  to  be 

the  larger  of  unity  or  1/ tan  2Sn.  Assuming  that  most  of  the  scattered  power 

from  a  single  scattering  event  is  lost  to  the  waveguide,  the  total  power  lost 
in  traversing  a  path  length  L  is  AP  ■  ELWId  *  ELdP.  The  waveguide  attenuation 
coefficient  is  therefore  AP/PL  ■  a  ■  Ed.  We  may  write 

£<J(0)  -  (l/2ir)  (2An/n)1/2  (af/6n2) .  (49) 
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Figure  9  shows  IOSA  dynamic  range  calculated  from  Eqs.  (48)  and  (25) 
as  a  function  of  a  with  <5n  as  a  parameter.  We  have  used  f  *  y  and  An  *  0.01. 

Values  of  n  less  than  0.05  are  not  plotted  because  for  these  values  the  scat¬ 
tered  light  cone  falls  within  the  acceptance  angle  of  the  waveguide,  and  the 
assumption  that  most  of  the  scattered  light  is  lost  to  the  waveguide  may  not 
be  valid.  Nevertheless,  for  plotted  values  of  <$n,  we  find  that  Mie  scat¬ 
tering  can  be  very  deleterious  to  IOSA  dynamic  range.  The  criterion  is  that 
the  scattering  centers  be  present  in  sufficient  numbers  to  produce  attenua¬ 
tion  levels  on  the  order  of  0.1  dB/cm. 
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Dynamic  Range,  dB 


Figure  9.  IOSA  dynamic  range  versus  waveguide  attenuation 
coefficient,  for  the  case  of  Mie  scattering 
centers. 
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III.  CONSIDERATION  OF  VARIOUS  SCATTERING  MECHANISMS  AND 
THEIR  INFLUENCE  ON  IOSA  DYNAMIC  RANGE 

The  results  of  the  previous  section  show  the  relative  dominance 
of  Rayleigh-Gans-Debye  and  Mie  scattering  to  Rayleigh  scattering  in  limiting 
IOSA  dynamic  range.  This  conclusion  is  supported  experimentally  by  results 
presented  in  the  following  section,  which  indicate  the  presence  of  small- 
angle  scattering  confined  to  within  a  few  degrees  of  the  beam  propagation 
direction. 

There  are  a  number  of  steps  in  the  fabrication  of  Ti-diffused 
LiNbO^  waveguides  where  either  Mie  or  RGD  scattering  centers  can  be  intro¬ 
duced.  The  purpose  of  this  section  is  to  identify  such  steps  and  evaluate 
whether  reasonable  estimates  of  scattering  center  size,  refractive  index, 
and  density  lead  to  unacceptable  values  of  IOSA  dynamic  range.  Candidate 
scattering  mechanisms  to  be  considered  are  (1)  surface  roughness,  (2)  lithium- 
titanate  surface  compounds,  (3)  the  occurrence  of  the  separated  phase  LiNb^Og, 
(4)  Ti-concentration  Inhomogeneity,  and  (5)  polishing  imperfections. 

SURFACE  ROUGHNESS 

Surface  roughness  scattering  centers  are  characterized  by  values 

of  <5n  on  the  order  of  unity  when  the  waveguide  supe  .ate  is  air.  A  more 

2 

careful  analysis  of  surface  roughness  scattering  shows  that  5n  ■  (n  -l)/2n  * 
0.87  for  Ti-diffused  LiNbO^  TE  modes.  We  will  also  identify  zQ  with  o,  the 
standard  notation  for  rms  surface  roughness,  and  a  with  the  roughness  auto¬ 
correlation  length.  We  will  continue  to  write  z  instead  of  o  to  avoid 

o 

confusion  with  our  notation  for  scattering  cross  section. 


Substituting  <Sn  ■  0.87  into  Eq.  (43),  we  find  the  relation  between 

z  and  <5n  to  be 
o  g 

-4  -1 

<5ng  ■  (5.59  x  10  pm  )zq  (50) 

for  a  step-index  waveguide.  In  the  case  of  surface  roughness  scattering  it 
is  possible  to  be  avoid  specifying  the  nature  of  the  waveguide  by  using 
the  differential  relationship  6ng  *  (3n  / 3D)zq.  For  LiNbO^  waveguides,  the 
derivative  is  accurately  described  by  the  relationship 

1/2 

3ng/3D  *  kQAn(2nAn)  '  (3b/3V),  (51) 

where  b  is  the  normalized  mode  index  and  V  is  the  normalized  depth.  The 
function  b(V)  is  described  by  the  universal  curve  appropriate  for  the  wave¬ 
guide  index  profile. In  the  case  of  a  Gaussian  profile,  3b/3V  *  0.1  for 
a  well  confined  mode  having  b  =  0.3  to  0.4.  The  use  of  this  value  in  Eq.  (51) , 
along  with  ■  0.83  pm,  n  =  2.2  and  An  ■  0.01,  results  in 

6n  «  (1.59  x  10  ^  pm~^)z  (52) 

g  o 

in  place  of  Eq.  (50)  . 

Figure  10  shows  the  corresponding  I0SA  dynamic  range  as  a  function 
of  zq  with  d  as  a  parameter.  Dynamic  range  levels  of  40  dB  are  achieved  for 
zq  10  A,  and  possibly  for  larger  values  of  zq  if  cl  is  as  small  as  several 
tens  of  microns.  National  Bureau  of  Standards  measurements  on  an  undiffused, 
polished  crystal  of  LiNbO^  obtained  from  Crystal  Technology  Corporation  show 
roughness  levels  _<  3.5  A  and  autocorrelation  lengths  in  the  range  20-40  pm.  We 
have  noted,  however,  using  differential-interference -phase-contrast  micros¬ 
copy,  an  increased  roughness  in  heavily  diffused  samples.  In  other  experi¬ 
ments,  we  have  measured  the  scattered  component  of  laser  light  reflected 

from  the  waveguide  surface,  and  used  the  data  to  evaluate  roughness  levels 

°  (11) 
of  about  30  A  and  autocorrelation  lengths  in  the  range  50-80  pm.  Accord¬ 
ing  to  Fig.  10  these  waveguides  would  be  associated  with  I0SA  dynamic  range 
values  of  about  32  dB,  or  8  dB  below  requirements. 

These  results  show  that  surface  roughness  can  be  a  factor  limiting 
IOSA  dynamic  range  to  substandard  values.  However,  the  potential  exists  for 
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repolishing  the  fabricated  waveguides  to  improve  dynamic  range.  Reduction 
of  zq  from  30  A  to  10  A  would  appear  to  be  sufficient.  We  have  measured  im¬ 
proved  scattering  performance  in  waveguides  that  were  polished  following  Ti 
diffusion,  as  indicated  in  Fig.  11.  However,  the  initial  waveguide  was 
heavily  diffused  with  a  720  A  thick  Ti  film.  When  waveguides  formed  with  Ti 
films  in  the  vicinity  of  200  A  were  polished  after  diffusion,  no  noticeable 
improvement  occurred.  The  efficacy  of  post-diffusion  polishing  for  waveguides 
of  the  type  fabricated  for  IOSA  applications  remains  to  be  certified. 

Li-Ti-0  SURFACE  COMPOUNDS 


Burns,  West,  Klein  and  Plew  have  presented  Secondary  Ion  Mass  Spec¬ 
troscopy  data  showing  the  presence  of  anomalously  high  Ti  concentrations 

(2) 

peaking  about  0.1  pm  beneath  the  surface  of  a  Ti-dif fused  LiNbO^  waveguide. 
This  peak,  shown  in  Fig.  12,  apparently  served  as  the  Ti  source  for  a 
Gaussian  profile  offset  from  the  surface  by  0.9  pm  and  having  a  diffusion 
depth  of  2  pm  (2.9  pm  relative  to  the  surface).  Holman  has  noted  that  Ti- 

diffused  LiNbO-  waveguides  will  occasionally  support  a  relatively  dim  mode 

J  (12) 
having  a  higher  effective  index  than  the  apparent  TEq  mode.  '  Bums  has 

suggested  that  this  mode  could  be  identified  with  the  near-surface  Ti  con- 
(13), 

centration  peak;  that  is,  the  energy  of  the  mode  could  be  localized 
largely  within  the  peak. 

We  observed  the  mode  in  the  sample  used  to  obtain  the  results  of 
Fig.  11.  At  the  time,  the  waveguide  had  been  polished  for  a  total  of  seven 

O 

hours  after  diffusion  of  a  720-A  Ti  film  for  3h  at  950°C.  The  mode  was  ob¬ 
served  for  TE  propagation  parallel  to  the  optical  axis  only.  The  effective 
index  of  the  mode  was  about  0.002  above  that  of  the  nominal  TE  mode,  and 
the  relative  scattered  intensity  at  a  given  angle  was  about  5  dB  greater  for 
the  dim  mode  than  for  the  nominal  mode.  After  11  hours  total  polishing,  the 
dim  mode  disappeared,  and  a  dramatic  improvement  in  waveguide  quality,  as 
shown  by  the  last  datum  in  Fig.  11,  ensued.  These  observations  suggest  that 
the  formation  of  near-surface  inhomogeneities  can  be  a  significant  source  of 
ln-plane  scattering,  apart  from  any  surface  roughness  generated  during  dif¬ 
fusion.  Bums,  Klein,  West  and  Flew  have  characterized  these  inhomogeneities 
as  Li-Ti-0  compounds  (U^TiO^  or  LijTi^O^)  formed  in  a  dilute  mixture  with 
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Figure  12.  Ti-concentration  versus  depth  beneath  the 
waveguide  surface.  From  Ref.  (2). 
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LiNbO^.  They  estimate  the  concentration  of  Li-Ti-0  compounds  to  be  less 

than  10  mol%. 

In  order  to  estimate  IOSA  dynamic  range  limitations  associated  with 
the  compounds,  it  is  necessary  to  specify  their  index  perturbation  Sn  and  their 
depth  beneath  the  surface  zq.  The  latter  quantity  is  estimated  from  Fig.  12 
as  0.4  pm.  The  former  quantity  depends  on  how  inhomogeneous ly  the  compounds 
are  distributed  within  zq.  The  index  perturbation  could  be  negligible  for  a 
uniform  distribution  or  it  could  be  as  large  as  the  index  difference  between 
the  Li-Ti-0  compounds  and  the  host  crystal.  The  index  of  L^TiO^  is  2.09,  and 
the  other  compounds  probably  have  similar  values,  suggesting  that  <5n  -  0.1  is 
a  good  estimate.  Since  the  compounds  are  believed  to  be  in  dilute  solution  of 
10  mol%  or  less  with  the  host  LiNbO^,  the  effective  value  of  6n  may  be  reduced 
another  order  of  magnitude  to  0.01. 

From  Eq.  (43),  based  on  an  analysis  for  a  step  index  waveguide,  the 

mode  index  perturbation  is  calculated  to  be  Sn  =  2.57  x  10  A  factor-of- 

three  increase  In  5n  is  probably  associated  with  a  Gaussian  waveguide  pro- 

© 

file,  as  we  found  to  be  the  case  for  the  surface  roughness  perturbation.  In 

this  event.  Fig.  7  shows  that  IOSA  dynamic  range  values  greater  than  40  dB 

require  a  3  pm.  This  indicates  reasonable  cause  for  concern  that  Li-Ti-0 

compound  formation  can  be  a  severe  problem  in  achieving  IOSA  dynamic  range 

specifications.  As  with  surface  roughness,  c.ie  remedy  Is  post  diffusion 

polishing  of  the  waveguide  to  remove  the  compound  layer.  It  is  also  possible 

that  the  layer  becomes  more  homogeneous  as  the  indiffusion  process  is  con- 

(14) 

tinued  for  longer  periods.  This  would  avoid  the  need  for  post-diffusion 

polishing  but  It  may  produce  weakly  guiding  or  multimode  waveguides  unsuited 
for  the  IOSA  application. 

Recently  we  have  noted  less  "orange  peel"  structure  on  the  surface 
of  waveguides  formed  in  the  presence  of  moist,  flowing  1^0  vapor. This 
orange  peel  is  probably  associated  with  the  increased  roughness  of  diffused 
waveguides  compared  to  the  initial  substrate.  It  may  also  be  indicative  of 
the  formation  of  subsurface  Li-Ti-0  compounds.  The  reduction  of  orange  peel 
using  moist  flowing  1^0  vapor  in  the  presence  of  oxygen  therefore  offers 
promise  for  reduced  waveguide  scattering. 


PHASE  SEPARATION 


Another  method  by  which  index  inhomogeneities  can  arise  in  LiNbO^ 
is  formation  of  the  separated  phase  LiNb^Og.  The  process  is  known  to  occur 
most  rapidly  when  the  crystal  is  heated  to  the  vicinity  of  800-850°C. 

Figure  13  shows  the  severe  surface  degradation  that  resulted  from  an  experi¬ 
ment  in  which  a  polished  sample  of  LiNbOg  was  held  at  850°C  for  10  minutes. 

The  phenomenon  of  phase  separation  may  have  served  to  decorate  strains  that 
were  generated  in  the  crystal  during  its  previous  history.  If  this  is  true, 
phase  separation  is  a  mechanism  for  enhancing  the  scattering  stength  of 
existing,  small  index  inhomogeneities . 

Although  heating  of  waveguide  samples  in  the  vicinity  of  850°C 
for  long  periods  does  not  normally  occur  during  conventional  fabrication 
procedures,  we  have  noticed  increased  scattering  from  some  waveguides  al¬ 
lowed  to  cool  slowly  from  the  diffusion  temperature.  This  prompted  us  to 
initiate  a  policy  of  quenching  waveguides  from  the  diffusion  temperature, 
usually  950°C,  to  600°C,  followed  by  a  slower  cool  to  room  temperature. 

This  resulted  in  waveguides  that  show  none  of  the  effects  of  Fig.  13.  How¬ 
ever,  crystals  from  a  subsequent  batch  fractured  during  this  quenching  opera¬ 
tion.  When  it  was  temporarily  abandoned  in  response,  there  was  no  apparent 
sacrifice  in  surface  quality.  We  conclude  that  problems  associated  with 
phase  separation  are  strongly  sample  dependent. 

Fabrication  techniques  that  minimize  induced  crystalline  strain 
may  be  the  best  guard  against  phase  separation  effects  in  LiNbOg  waveguides. 
These  would  probably  include  techniques  which  reduce  Li-Ti-0  compound  inhomo¬ 
geneities  near  the  waveguide  surface  as  well.  To  the  extent  that  phase  sepa¬ 
ration  occurs  near  the  vicinity  of  these  inhomogeneities,  formation  of  LiNb^Og 
may  add  to  the  inherent  scattering.  The  refractive  indices  of  the  monoclinic 
crystal  LiNb^Og  are  2.28,  2.36,  and  2.40.  The  scattering  effect  from  this 
phase  is  therefore  comparable  to  that  from  similarly  sized  lithium-titanate 
inhomogeneities . 

In  event  that  LiNb^Og  aggregates  extend  through  the  waveguide  layer 
into  the  substrate,  they  will  likely  produce  either  Mie  scattering  centers 
with  5n<  0.1  or  RGB  scattering  centers  with  6n  >>10~^.  In  either  situation 

-v  g 

IOSA  dynamic  range  values  less  than  the  required  40  dB  should  result.  In  the 
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L3.  Nomarski  micrograph  showing  the  results  of 
a  10  rain  heat  treatment  of  LiNb03  at  850°C 
in  flowing  Oj.  (500X) 
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worst  case,  phase  separation  produces  surface  degradation  like  that  shown  in 
Fig.  13.  The  index  inhomogeneity  associated  with  these  scattering  centers 
6n  -  1.  If  their  depth  is  0.4  ym,  they  will  produce  40  dB  more  scattered 
intensity  than  we  estimate  for  near  surface  Li-Ti-0  compounds.  While  this 
would  be  disastrous  for  the  IOSA  application,  we  note  that  surface  degrada¬ 
tion  like  that  of  Fig.  13  is  not  the  norm  for  TitLiNbO^  waveguides,  and  it 
can  be  readily  seen  when  it  occurs. 


Ti  CONCENTRATION  INHOMOGENEITIES 


One  possible  source  of  scattering  in  Ti:LiNbO^  waveguides  is  non- 
uniform  Ti  concentration  resulting  from  defects  in  the  initial  Ti  film.  For 
example,  dust  trapped  beneath  the  film  during  its  formation  can  prevent  dif¬ 
fusion  of  the  metal  into  the  waveguide  surface.  This  could  be  caused  by  an 
explosive  removal  of  the  film  in  the  vicinity  of  the  dust  particle  as  the 
high  diffusion  temperature  is  approached.  Alternatively,  the  dust  could 
carbonize,  react  with  the  metal  or  the  crystal,  or  otherwise  inhibit  the 
normal  diffusion  mechanism.  Figure  14  shows  a  micrograph  of  the  edge  of  a 
nominally  good  waveguide  where  both  of  these  mechanisms  may  be  in  evidence. 

The  smooth  round  artifacts  that  have  the  appearance  of  craters  are  consistent 
with  the  explosive  removal  hypothesis,  in  that  Ti  is  known  to  produce  a  thick¬ 
ness  change  in  diffused  LiNbO^.  Where  this  change  is  not  observed,  we  may 
not  have  diffused  Ti.  In  contrast  globular  aggregates  in  Fig.  14  may  show 
regions  where  chemical  reactions  have  inhibited  the  diffusion  mechanism. 

None  of  the  structures  in  Fig.  14  are  observed  in  the  center  of  the 
waveguide  sample,  where  surface  cleanliness  is  easier  to  achieve  and  maintain. 
Nevertheless  an  occasional  dust  particle  could  produce  a  Ti-free  region  of 
the  waveguide  where  the  effective  index  is  reduced  from  ng  to  the  substrate 
index  n.  For  a  well-confined  TEq  mode  in  a  single-mode  waveguide,  the  re¬ 
sultant  index  inhomogeneity  would  be  6n  2  0.003.  While  this  is  a  fairly 
strong  RGD  scattering  center  or  a  weak  Mie  scattering  center,  IOSA  dynamic 
range  will  only  be  seriously  impaired  if  the  density  of  these  scattering 
centers  is  sufficiently  large.  Figure  15  shows  a  plot  of  IOSA  dynamic 
range  versus  Z  with  scattering  center  size  a  as  a  parameter.  The  plot  is 
based  on  an  RGD  scattering  analysis.  The  size  of  the  smooth  craters  in 
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Fig.  14.  Nomarski  micrograph  showing  artifacts 
near  the  edge  of  a  Ti-diffused  LiNbC>3 
waveguide.  (500X) 
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Fig.  14  is  in  the  range  3-10  um.  Therefore,  achievement  of  an  IOSA  dynamic 

range  greater  than  40  dB  requires  that  the  density  of  such  craters  be  reduced 

-2  -2 

to  about  0.1  mm  or  10  cm  .  This  is  within  the  realm  of  careful  sample 
preparation  procedures. 


POLISHING  IMPERFECTIONS 


When  glass  is  polished  improperly,  it  is  possible  to  obtain  a  smooth 
surface  on  top  of  a  layer  containing  voids  and  inhomogeneities.  This  occurs 
because  of  the  ability  of  glass  to  flow  during  the  polishing  operation,  and 
thereby  to  cover  as  well  as  remove  damage  generated  during  previous  polishing 
steps.  Although  the  polishing  characteristics  of  LiNbO^  are  very  different 
from  glass,  there  is  concern  that  a  similar  situation  could  exist  in  the  case 
of  Improperly  polished  LiNbO^- 

Subsurface  damage,  if  it  exists,  must  be  comprised  of  scattering 
centers  small  compared  to  a  wavelength  or  it  would  very  likely  be  visible 
under  a  microscope.  If  the  damage  is  small  compared  to  a  wavelength,  IOSA 
dynamic  range  is  described  by  Fig.  5  for  Rayleigh  scattering.  The  achieve¬ 
ment  of  40  dB  or  more  dynamic  range  is  tied  to  the  reduction  of  waveguide 
attenuation  levels  to  0.25  dB/cm  or  less.  Since  this  is  within  the  state 
of  the  art,  and  since  our  experiments  to  be  described  do  not  show  the  pre¬ 
sence  of  Rayleigh  Scattering,  we  do  not  perceive  polishing  imperfections  to 
be  a  serious  problem  limiting  IOSA  dynamic  range. 

SUMMARY 


The  analysis  of  this  section  shows  that  RGD  scattering  centers 
associated  with  surface  roughness  and  near-surface  compound  formation  should 
be  among  the  hardest  to  diagnose  and  prevent,  as  well  as  among  the  most 
damaging  to  IOSA  dynamic  range.  Much  of  our  experimental  evidence  in  the 
following  section  will  be  found  to  support  this  assertion. 
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IV.  EXPERIMENTAL  STUDIES  OF  IN-PLANE 
SCATTERING  IN  TirLiNbOj  WAVEGUIDES 

MEASUREMENT  OBJECTIVE 


The  objective  of  our  experimental  program  has  been  to  measure 
the  parameters  Io(0)  and  q  using  observations  of  in-plane  scattering  in 
diffused  Ti:LiNbO^  waveguides  fabricated  in  various  ways.  According  to 
the  theory  developed  in  Sec.  II,  measurement  of  Io(0)  and  q  allows  us  to 
predict  the  scattering  limits  to  dynamic  range  in  the  IOSA  application. 

By  correlating  the  values  of  the  scattering  parameters  with  waveguide 
fabrication  parameters,  we  may,  in  addition,  hope  to  develop  a  set  of 
optimized  fabrication  conditions.  We  may  also,  by  careful  interpretation 
and  analysis  of  results ,  hope  to  reach  valid  conclusions  about  the  dominant 
source  of  scattering  in  the  waveguides. 

The  obvious  method  for  measuring  the  in-plane  scattering  parameters 
is  shown  in  Fig.  16.  A  collimated  laser  beam  is  coupled  into  and  out  of  the 
waveguide  and  brought  to  a  focus.  In-plane  scattered  energy  is  distributed 
along  the  "m-line"  consisting  of  a  brightly  focused  central  spot  and  a  dif¬ 
fuse  component  that  decreases  in  intensity  symmetrically  with  distance  on 
either  side.  This  component  fluctuates  in  intensity  owing  to  optical  inter¬ 
ference  of  coherent  rays  generated  by  randomly  positioned  scattering 
centers. 

If  a  slit  of  width  s  is  positioned  a  distance  y  from  the  focus, 
the  slit  will  pass  rays  scattered  in  the  waveguide  plane  within  an  angular 
sector  extending  from  to  $  4-  A$/2,  where 

<j>  -  y/nf, 

i 

*  s/nf  .  (53 

The  focal  length  of  the  lens  is  f ,  and  a  factor  n  is  introduced  to  account 
for  the  refraction  of  light  rays  in  leaving  the  waveguide  medium  and  enter¬ 
ing  air.  From  Eq.  (2),  the  scattered  power  passed  by  the  slit  is 


AP(y)  -  PLlo(y/nf) (s/nf) . 


(54 


At  y  ",  0*  the  scattered  power  is  swamped  by  the  focused  power  associated 
with  the  unscattered  beam.  This  is 

P(0)  -  P  erf (tt21^2  sW/4f A  ) .  (55) 

In  processing  the  experimental  data,  the  scattered  power  is  usually  divided 
by  the  peak  power.  The  measured  quantity  is  then 

AP(y)/P(0)  -  LZa(y/nf)(s/nf)  [erf(Tr21/2  sW/4fA  i]-1.  (56) 

For  a  slit  width  small  compared  to  the  focal  spot  size,  this  result  reduces  to 

1/2 

AP($)/P(0)  =  (2/ tt) X/  lL/W)  (AQ/n)Zo  (<j>) ,  (57) 

where  we  have  reinserted  $  in  place  of  y/nf.  All  parameters  in  the  equation 
other  than  So (<p)  are  known.  Therefore  the  data  is  directly  proportional  to 
the  scattering  parameters  and  the  constant  of  proportionality  can  be  calculated. 

THE  DEVELOPMENT  OF 
EXPERIMENTAL  PROCEDURES 


Evidence  for  Scattering  from  Prism  Coupling 

The  main  experimental  problem  is  that  of  insuring  that  the  measured 
scattered  light  is  produced  only  by  waveguide  scattering  along  the  path  L  be¬ 
tween  input  and  output  coupling.  All  optical  components  that  refract,  re¬ 
flect,  or  aperture  the  beam  contribute  to  light  intensity  outside  the  focal 
spot.  This  is  effectively  a  source  of  measurement  noise  that  must  be  either 
eliminated  or  carefully  accounted  for.  Most  optical  components  can  be 
selected  to  be  of  sufficiently  high  optical  quality  that  scattering  from 
them  is  not  a  serious  problem.  In  addition,  scattered  light  from  components 
well  removed  from  the  measurement  plane  can  be  effectively  shielded  from  the 
plane  by  careful  aperturing. 

The  source  of  the  worst  problems  is  scattering  generated  during 
input  and  output  coupling  from  the  waveguide.  This  source  of  scattering 
is  in-plane  and  cannot  effectively  be  apertured  because  of  its  close  proxim¬ 
ity  to  the  scattered  light  from  the  waveguide.  In  the  case  of  prism  coupl¬ 
ing,  evidence  for  coupling-related  scattering  is  the  observation  of  nearly 
identical  m-lines  for  different  coupling-prism  separations.  Figure  17  shows 
one  example  of  the  variation  of  total  scattered  power  as  a  function  of 
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Figure  17.  Relative  scattered  power  vs.  coupling  prism  separa¬ 
tion  for  Ti:LiNb03  waveguide  //166  (TE  mode,  k|jc). 


coupling  prism  separation.  The  experiment  measured  the  total  power  scat¬ 
tered  into  an  angular  sector  about  1.3°  on  either  side  of  the  optical  axis. 

The  results  were  normalized  by  dividing  by  the  unscattered  power  leaving 
the  waveguide.  For  RGD  scattering,  the  expected  variation  is  obtained  by 
inserting  Eq.  (1)  into  Eq.  (2)  and  integrating  from  $  ■  -1.3°  to  $  *  +1.3°. 

The  result  is 

AP  /P  "  2LZo(0J (koa)_1  tan~1(0.024  kQa) ,  (58) 

A  plot  of  AP  /P  versus  path  length  L  should  plot  as  a  straight 
line  with  zero  intercept.  The  results  of  Fig.  17  show,  in  contrast,  a 
finite  intercept  associated  with  prism-related  scattering.  The  best  fit¬ 
ting  straight  line  Indicates  that  at  a  path  length  of  15  mm,  only  9%  of 
the  total  scattering  observed  comes  from  the  waveguide.  91%  comes  from 
prism  coupling.  The  dashed  straight  lines  are  within  one  standard  devia¬ 
tion  of  the  best  fitting  straight  line,  and  they  show  the  possibility  that 
the  correct  straight  line  through  the  data  has  zero  slope,  indicative  of 
the  result  £o(0)  »  0.  That  is,  the  waveguide  is  scatter  free. 

The  best  straight  line  fit  is  used  to  obtain  the  lower  curve  in 
the  graph  of  Fig.  18.  This  curve  shows  that  the  waveguide  used  to  obtain 
the  Fig.  17  results  is  capable  of  providing  an  IOSA  dynamic  range  better 
than  40  dB  if  scattering  center  size  is  greater  than  about  3  pm.  However, 
the  upper  curve  in  Fig.  18  is  associated  with  the  dashed  fit  to  the  data  of 
Fig.  17  showing  a  slope  just  one  standard  deviation  higher  than  the  best 
fitting  slope.  And  this  curve  indicates  that  a  dynamic  range  of  40  dB  is 
not  achieved  for  any  scattering-center  size  in  the  range  from  1-100  pm. 

Data  taken  using  the  TM  mode  instead  of  the  TE  mode  are  also 

o  o 

disheartening.  For  that  data,  the  most  likely  IOSA  dynamic  range  curve  is  the 
same  as  the  top  curve  in  Fig.  18.  A  further  3  dB  dynamic  range  is  sacrificed 
when  the  steepest  slope  within  one  standard  deviation  of  the  best  slope  is 
employed  in  the  calculation.  However  the  shallowest  slope  within  one  stan¬ 
dard  deviation  of  the  best  slope  allows  for  the  possibility  of  zero  wave¬ 
guide  scattering.  In  short,  the  data  tells  us  that  the  waveguide  may  be 
adequate  for  the  IOSA  application,  or  that  10  dB  further  scattering  reduc¬ 
tion  is  needed.  Unfortunately  the  error  bars  do  not  allow  for  a  more  pre¬ 
cise  conclusion. 


Dynamic  Range,  dB 


Figure  18.  IOSA  dynamic  range  versus  scattering  center  size, 
based  on  the  experimental  data  of  Figure  17. 


54 


Figure  19  shows  a  different  type  of  measurement  on  the  same  wave¬ 
guide,  in  which  10  log  AP(<f>)/P(Q)  Is  plotted  as  a  function  of  scattering 
angle  in  the  waveguide.  The  measurement  was  made  by  continuously  scanning 
a  slit  along  the  focused  m-line.  A  fiber  bundle  coupled  to  the  slit  col¬ 
lected  the  light  passed  by  it  and  transmitted  the  light  a  photomultiplier 
tube.  Figure  19  actually  contains  the  results  of  two  separate  experiments, 
in  which  the  coupling-  prism  separation  was  either  15  mm  or  1  mm,  approxi¬ 
mately.  The  data  for  the  two  experiments  are  superimposed,  with  the  dashed 
curve  representing  the  results  obtained  using  closely  spaced  coupling  prisms. 
The  overall  envelopes  of  the  two  curves  are  nearly  identical.  A  statistical 
analysis  has  been  used  to  generate  Lorentzian  shaped  best  fitting  curves  for 
the  two  measurements.  These  are  also  plotted  in  Fig.  19.  They  indicate  a 
slight  excess  of  scattering  for  the  case  of  well-separated  coupling  prisms. 
However  the  excess  is  small  compared  with  the  fluctuations  associated  with 
laser  speckle.  Accordingly  it  is  difficult  to  use  Eq.  (55)  to  obtain  an 
accurate  determination  of  £a(4>),  as  required  to  predict  IOSA  dynamic  range. 

To  the  extent  that  the  data  can  be  trusted,  the  two  Lorentzian  curves  deter- 

-2  -1 

mine  the  values  Ecr(0)  ■  5.7  x  10  mm  and  a  *  20  pm.  The  corresponding 
IOSA  dynamic  range  at  the  fourth  off-axis  pixel  is  32  dB.  This  is  h  dB 
worse  than  the  more  pessimistic  value  indicated  in  Fig.  18. 

Mechanisms  for  Scattering  from  Prism  Coupling 

We  have  considered  two  mechanisms  by  which  prism  coupling  can  add 
to  the  amount  of  in-plane  scattered  light,  even  if  the  prism  contains  no 
scattering  centers  itself. 

One  is  based  on  the  possibility  that  the  electromagnetic  field 
at  the  surface  is  enhanced  in  the  region  where  light  is  being  prism 
coupled  into  the  waveguide.  The  field  in  this  region  consists  of  the  wave¬ 
guide  modal  field  plus  the  evanescent  field  generated  in  the  base  of  the 
coupling  prism  by  the  incoming  laser  beam.  Some  of  this  field  is  reflected 
and  does  not  contribute  to  the  modal  field.  Yet  it  is  phase-matched  with 
the  modal  field  and  can  presumably  generated  scattered  light  which  is 
trapped  in  the  waveguide.  This  scattering  mechanism  relies  on  the  incoming 
laser  beam  and  so  it  exists  only  for  prism  input  coupling.  We  have  observed 
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Figure  19.  In-plane  scattered  energy  distribution  for  sample  166 
(TM  ,  kjc) •  Scattering  angles  are  referenced  to  wave¬ 
guide  material. 


cases  for  which  light  was  end-fire  coupled  into  a  waveguide  and  prism  coupled 
out.  Qualitatively,  the  scattering  in  this  situation  was  less  than  when  the 
coupling  roles  were  reversed.  However,  data  to  support  this  were  not  acquired. 
Calculations  presented  in  Ref.  11  show  that  the  postulated  mechanism  is 
realistic,  but  they  do  not  support  the  extent  to  which  prism  scattering 
appears  to  dominate  waveguide  scattering  in  practice. 

A  second  mechanism  for  prism  scattering  is  found  in  spatial 
variations  in  coupling  efficiency  across  the  width  of  an  incoming  beam.  After 
coupling,  the  profile  of  an  initially  Gaussian  beam  may  contain  structure  of 
characteristic  size  A.  The  diffraction  pattern  produced  in  the  far  field 
is  indistinguishable  from  that  produced  by  scattering  centers  of  the  same 
size.  This  mechanism  exists  for  both  input  and  output  prism  coupling. 

To  estimate  the  size  of  this  effect,  let  us  assume  that  A(y)  is 
the  amplitude  profile  of  an  ideal  Gaussian  beam,  y  being  the  transverse 
direction  in  the  plane  of  the  waveguide.  The  actual  profile  of  an  in- 
coupled  beam  is  [1  +  t(y)]A(y),  where  t(y)  is  related  to  the  spatial 
variation  of  input  coupling  efficiency.  If  the  beam  of  light  is  coupled 
out  of  the  waveguide  and  focused  by  a  lens  of  focal  length  f,  the  field 
at  the  focal  plane  will  be 


E(Y)  -  Eo(Y)  +  E  (Y) 

00 

E  (Y)  -  (ifAor1/2  /  t(y)A(y)e'27riyY/fA°  dy  (59) 

^  —00 

where  E  (Y)  is  the  ideal  Gaussian  field  and  E  (Y)  is  the  field  associated  with 

o  p 

nonuniform  prism  coupling  efficiency.  The  relative  power  in  the  field  Ep 

is  AP  (Y)/P(0)  *  <E  E  >/E2(0) .  The  bracketed  quantity  is 
P  p  p  o 

<E  E  >  -  (fXQ)  1  /  /  dydy’  <t(y)t(y')>  A(y)A(y')e-2lTi^y  y  ^Y  fAo  (60) 

A  standard  approach  to  integrals  of  this  type  is  to  set  <t(y)t(y')>  ■ 

2 

tQ  exp(-|y-y' |/A) ,  where  t  is  the  rms  amplitude  transmission  fluctuation, 
and  A  is  the  correlation  length  associated  with  these  fluctuations.  If  the 
beam  waist  W  is  large  compared  to  A,  we  find 
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(61) 


APp(4>)/P(0) 


2<2/ir) 


1/2  2 
t 

o 


(A/W)  [1  +  (kQA4>)  ] 


where  Y/f.  The  angular  variation  Is  a  Lorentzlan  function  In  which  A 
plays  the  same  role  chat  the  scattering-center  size  a  plays  in  scattering 
analysis . 

To  estimate  the  magnitude  of  this  pseudo-scattering  effect,  we  have 
measured  output-coupled  beam  profiles  within  2  cm  of  the  output  coupling 
prism.  The  result  is  shown  in  Fig.  20.  The  detector  was  a  linear  photo¬ 
diode  array  with  a  25-um  center- to-center  pixel  spacing.  The  deviations 

of  the  beam  profile  from  an  ideal  Gaussian  shape  are  described  by 
2 

[1  +  t(y)]  -  1  +  2t(y).  Without  being  quantitative,  an  examination  of 

Fig.  20  suggests  that  fluctuations  from  a  smooth,  average  beam  profile  are 
typically  30%,  and  that  individual  fluctuations  span  about  4  channels  or 

100  pm.  Thus  we  estimate  2t  -  0.30  and  A  *  100  pm.  The  on-axis  scattered 

°  1/2 
power  associated  with  prism  coupling  effects  is  therefore  2(2/ir)  t  (A/W) 

-4  ° 

*  3.6  x  10  or  -24.4  dB  when  W  *  1  mm.  This  is  comparable  to  what  we  have 

observed  in  many  experiments.  The  on-axis  scattered  power  from  the  wave¬ 
guide  is  (L/W)lcr(0) XQ/n.  If  this  value  is  small  compared  to  -25  dB,  the  wave¬ 
guide  scattering  is  swamped  by  the  pseudo-scattering  associated  with  input 
coupling  variations. 


Experimental  Approaches  to  the  Prism-Scattering  Problem 

To  circumvent  this  problem  we  have  utilized  five  different  experi¬ 
mental  approaches:  (1)  Acquire  scattering  data  using  two  or  more  separations 
between  the  input  and  output  coupling  prisms.  (2)  Avoid  the  use  of  prisms 
by  implementing  end-fire  entry  and  exit  of  light  from  polished  waveguide 
terminations.  (3)  Use  waveguides  with  a  long  path  length  to  increase  the 
amount  of  waveguide  scattering  detected.  (4)  Take  data  at  increasingly 
large  scattering  angles  to  avoid  interference  with  prism  scattering  that 
predominates  at  small  angles.  (5)  Use  an  imaging  technique  to  avoid  detec¬ 
tion  of  light  scattered  during  input  coupling.  Our  comments  on  these 
approaches  ere  as  follows: 

(1)  Vary  coupling  prism  separation.  We  have  discussed  this  ap¬ 
proach  in  connection  with  Figs.  17-19  and  we  will  use  it  again  later  with 
reference  to  other  data. 
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Figure  20.  Output  beam  profile  showing  effects  of 
nonuniform  prism  coupling  efficiency. 
(Pixel  spacing  is  25  um) . 
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(2)  Use  end-fire  coupling.  With  this  method  we  have  found  scatter¬ 
ing  levels  greater  than  obtained  using  prism  coupling. 

(3)  Increase  waveguide  path  length.  We  have  used  paths  as  long  as 
57  mm  without  increasing  the  waveguide  scattering  level  significantly  above 
the  background  level. 

(4)  Examine  larger  scattering  angles.  If  A  is  larger  than  a,  then  the 

prism-scattered  component  described  by  Eq.  59  will  eventually  fall  below  the 

waveguide  scattered  component  described  by  Eq.  (55).  To  take  advantage  of 

this  we  have  increased  the  angular  range  of  our  measurements  a  factor  of  about 

five  beyond  the  range  indicated  in  Fig.  19.  At  large  angles,  however,  the 

2  2  2 

relative  scattered  power  varies  as  Za(0)/k  a  $  .  The  data  serves  only  to 

_2 

determine  the  combination  of  parameters  £o(0)<2  .  Not  enough  information  is 

available  to  determine  Eo(0)  and  CL  separately,  as  required  for  IOSA  dynamic- 
range  prediction. 

(5)  Use  of  an  imaging  technique.  In  this  technique  a  focused  beam  is 
coupled  into  the  waveguide  using  either  prism  or  end-fire  methods.  The  beam 
diverges  as  it  propagates  in  the  waveguide.  After  output  coupling  the  beam 
is  again  brought  to  a  focus.  Light  scattered  in  the  region  of  input  coupling 
is  imaged  at  the  focus.  Consequently  it  is  swamped  by  unscattered  light  and 
does  not  contribute  to  the  measured  level  of  scattered  power.  This  technic ue 
is  useful  only  if  output  coupling  can  be  performed  without  adding  scattered 
light  to  the  beam.  Thus,  it  is  not  useful  If  coupling  efficiency  nonuni¬ 
formities  are  the  dominant  cause  of  prism-associated  scattering.  Possibly 
for  this  reason,  our  results  with  this  technique  were  not  encouraging. 

EXPERIMENTAL  RESULTS 

Comparison  of  Waveguide  Fabrication  Methods 

The  samples  employed  in  this  research  were  of  two  general  categories : 
those  fabricated  at  Battelle  and  those  fabricated  at  Westinghouse.  In  both 
cases,  polished  crystals  of  LiNbO^  were  purchased  from  Crystal  Technology. 

The  Battelle  waveguides  were  formed  on  25  mm  x  25  mm  x  3  mm  substrates, 
while  the  Westinghouse  waveguides  were  formed  on  75  mm  x  25  mm  x  3  mm  sub¬ 
strates.  The  Westinghouse  waveguides  were  made  identical  to  those  used  in 

(4) 

the  fabrication  of  their  IOSA  systems.  Titanium  film  thicknesses  were  in  the 


range  275-300  A.  Diffusions  were  performed  at  1000°C  for  6  hr  in  an  environ¬ 
ment  containing  LiNbO^  powder,  to  prevent  outdiffusion,  and  sufficient  0^  to 
prevent  discoloration.  The  Battelle  samples  were  made  using  Ti-film  thick- 

O 

nesses  in  the  range  175-600  A,  but  most  samples  were  made  with  thicknesses 
about  equal  to  those  used  by  Westinghouse.  Battelle  diffusion  temperatures 
tended  to  be  similar  (<  1000°C)  and  shorter  (typically  3  hr.).  Contrary  to 
the  common  experience,  we  find  these  conditions  sufficient  to  complete  the 
diffusion  process.  This  possibly  has  to  do  with  the  initial  oxygen  content 
of  the  sputtered  or  e-beam  evaporated  Ti  film.  However,  one  would  antici¬ 
pate  that  the  Battelle  waveguides  are  shallower  and  have  a  larger  surface 
index  change.  According  to  the  discussion  of  Sec.  Ill,  this  should  produce 
greater  scattering  from  surface  roughness  or  surface  index  inhomogeneities. 
However,  we  have  observed  that  recent  Battelle  fabrication  methods  employing 
flowing  water  vapor  produce  less  surface  disruption  than  the  above- 

described  Westinghouse  procedure. 

Methods  of  Data  Acquisition 

Figure  21  shows  one  of  two  experimental  configurations  used  to  measure 
in-plane  scattering  parameters  in  Ti:LiNb0^  waveguides  provided  by  Westing¬ 
house.  The  system  uses  the  conventional  Fourier-transform  geometry  in  which 
a  collimated  guided  beam  is  focused  by  a  lens  to  provide  the  far  field  dis¬ 
tribution  from  which  scattering  parameters  are  obtained  in  accord  with  Eq.  (57). 
In  Fig.  21  a  second  lens  is  employed  to  produce  a  25X  magnified  image  of  the 
far-field  distribution  on  the  linear  detector  array  of  an  optical  mutichannel 
analyzer  (OMA) .  The  length  of  the  detector  array  is  12.7  mm.  Thus  the 
array  records  the  intensity  distribution  over  a  length  12.7  mm/25  =  0.5  mm 
measured  in  the  focal  plane  of  the  first  lens.  Since  the  focal  length  of 
the  first  lens  is  45  mm,  the  detector  array  subtends  an  angular  range  of 
0.5/45  or  0.64°.  In  order  to  examine  a  total  angular  range  of  6°,  it  is 
necessary  to  step  the  second  lens  perpendicular  to  the  beam  through  approxi¬ 
mately  10  positions  separated  by  0.5  mm.  At  each  position,  a  photograph  of 
the  OMA  output  constitutes  part  of  the  raw  data.  One  such  photograph  is 
shown  in  Fig.  22. 
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The  OMA  is  calibrated  so  that  the  scattered  power  into  any  of  the 

500  pixels  may  be  inferred  from  the  height  of  the  signal  measured  from  the 

photograph  or  read  from  the  instrument  display  panel  when  the  data  is  taken. 
Although  the  peak  power  from  the  focused  unscattered  beam  is  shielded  from 
the  detector  array,  this  power  is  known  from  a  separate  measurement  of  the 
total  power  in  the  beam,  plus  an  initial  determination  of  the  fraction  of 
this  power  captured  by  a  centrally  illuminated  pixel  when  the  beam  is  not 
shielded  from  the  array. 

Thus,  AP(4>)/P(0)  may  be  measured  for  scattering  angles  ranging  from 
4>  =  0°  to  6°  in  air,  or  0°  to  6° 12.2  =  2.7°  in  the  waveguide  medium.  Since 
the  scattered  power  goes  down  as  the  scattering  angle  <j>  increases,  the  laser 
power  is  increased  as  needed  to  insure  that  the  level  of  scattered  power  in¬ 
cident  on  the  detector  array  is  maintained  at  a  suitable  value  in  excess  of 

the  noise  level  but  below  saturation.  We  have  measured  and  plotted  the  scat¬ 
tered  power  only  in  the  fluctuation  peaks  in  order  to  provide  the  envelope 
of  the  actual  distribution.  For  data  like  that  shown  in  Fig.  22  we  select 
three  to  five  roughly  equally  spaced  fluctuation  peaks  for  use  in  the  analysis. 

Westinghouse  Samples 

Westinghouse  Corporation  provided  us  with  two  samples  made  by  similar 
techniques  but  of  very  different  quality.  Most  of  our  results  were  obtained 
with  the  better  of  the  two  samples,  called  Wl.  Using  the  experimental  arrange¬ 
ment  of  Fig.  21,  we  took  data  with  Wl  for  prism  separation  of  about  53  mm  and 
about  1  mm,  and  plotted  the  results  on  the  same  graph.  Several  months  later 
we  repeated  the  experiment  to  check  for  reproducibility.  Figures  23  and  24 
show  the  results  of  these  two  runs. 

Data  obtained  with  input  and  output  prisms  separated  by  only  1  mm 
are  plotted  as  plus  signs,  while  data  obtained  with  prisms  separated  by  52 
or  53  mm  are  plotted  as  dots.  Scattering  contributed  from  the  waveguide 
path  should  cause  the  dots  to  be  located  above  the  plusses,  assuming 
that  the  scattering  associated  with  prism  coupling  is  the  same  for  both 
prism  separations.  We  find  that  dots  lie  above  plusses  only  for  scattering 
angles  less  than  3°  in  Fig.  23.  In  Fig.  24,  the  scattered  power  at  similarly 
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Scattering  Angle  in  Air,  Degrees 

are  23.  Normalized  scattered  power  versus  scattering  angle  for  the  TMq  mode  of  Sample  Wl. 

Dots  correspond  to  53  mm  prism  separation.  Pluses  correspond  to  less  than  1  mm 
prism  separation. 


Scattering  Angle  in  Air,  Degrees 

Normalized  scattered  power  versus  scattering  angle  for  the  TMq 
mode  of  sample  Wl.  Dots  correspond  to  52  mm  prism  separation. 
Pluses  correspond  to  less  than  1  mm  prism  separation.  Data 
were  taken  two  months  after  data  of  Fig.  23  using  essentially 
the  same  experimental  set  up. 


small  angles  Is  nearly  Indpendent  of  prism  separation.  At  large  angles, 
plusses  are  above  the  dots  In  both  figures,  Indicating  the  unphysical  result 
that  scattering  is  greater  when  waveguide  path  length  Is  less. 

We  later  learned  that  background  scattering  from  optical  components 
other  than  prisms  was  dominant  at  these  large  scattering  angles.  This  scat¬ 
tered  component  was  apparently  detected  in  greater  amounts  when  the  prisms 
were  placed  close  together,  owing  to  the  experimental  geometry.  Because  of 
it,  we  have  to  ignore  the  data  at  scattering  angles  larger  than  2.2°  in 
both  Fig.  23  and  24.  At  very  small  scattering  angles,  less  than  0.5°,  we 
remain  unsure  of  the  data  because  of  the  pseudo-scattering  associated  with 
spatial  variations  in  coupling  efficiency,  discussed  earlier. 

As  a  result,  we  concentrate  our  analysis  on  data  falling  in  the  range 
0.5°-2.2°.  Moreover,  we  correct  the  data  in  this  range  by  subtracting  that 
part  of  the  scattered  power  contributed  by  the  nearly  constant  background 
at  larger  angles.  This  background  was  estimat ’d  at  -56  dB  for  the  +  data 
and  -58  dB  for  the  •  data.  The  results  are  found  to  closely  fit  a  Lorentzian 
angular  distribution  characteristic  of  RGD  scattering  centers, 

AP(*)/P(0)  -  (AP/P)  [1  +  (k  air/180)2^2]'1  (62) 

UlaA  O 

where  (AP/P)mny  depends  on  waveguide  and  prism  scattering  parameters.  Figures 

25  and  26  show  the  corrected  data  from  Figs.  23  and  24,  respectively,  plotted 

2 

in  the  form  P(0)/AP(<{i)  versus  $  .  When  plotted  in  this  form,  Eq.  (62)  becomes 

a  straight  line.  Fits  to  the  data  shown  in  the  figures  are  determined  using 

the  method  of  least  squares.  Table  IV  shows  the  tabulated  values  for  slope 

and  intercept  along  with  the  corresponding  standard  deviations. 

Ideally,  the  measured  slope  and  intercept  values  serve  to  determine 

(AP/P)^^  and  a.  However,  the  negative  value  for  the  intercept  found 

in  all  cases  is  unphysical.  The  experiment  simply  cannot  be  used  to  determine 

the  on-axis  scattering  to  any  useful  degree  of  accuracy.  The  best  that 

can  be  done  is  to  use  the  slope  information  to  obtain  the  combination  of  parameters 

2 

(AP/P)  /(k  air/180)  .  If  we  refer  to  this  combination  as  Q(L) ,  where  L 
max  o 

is  the  coupling  prism  separation,  then  the  waveguide  contribution  to  scatter¬ 
ing  at  large  angles  is 
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TABLE  IV.  SCATTERING  PARAMETERS  FOR  WESTINGHOUSE  WAVEGUIDE  W-l 
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In  this  expression  L  is  52  or  53  mm  and  L  .  is  1  mm.  Comparison  of 

max  min  _ 

Eq.  (63)  with  our  experimental  results  enables  us  to  calculate  Eo(0)/a  . 

2 

For  the  first  run  made  with  waveguide  W  1,  we  find  Ecr(0)/<z  *  6.18  x  10  pm 

2  —8 

For  the  second  run,  we  obtain  the  unphysical  result  Zo(0)/a  =  -1.04  x  10 

_3 

pm  .  However,  if  we  use  Q-values  removed  one  standard  deviation  from  the 

2  -8 

best  value,  we  can  obtain  a  positive  result;  namely  Za( 0)/a  =*  0.84  x  10 

-3  2 

pm  .  The  smallest  value  of  Io(0 )/a  obtained  for  the  first  run  using  Q-values 

2  —8 

within  one  standard  deviation  of  the  best  values  in  Za(0)/a  =  3.88  x  10 

_3 

pm  .  The  inconsistencies  in  the  two  runs  may  result  from  changes  in  the 
prism  contribution  to  scattering  resulting  from  different  prism  placements. 

Figure  27  shows  the  predicted  IOSA  dynamic-range  variation  with  scat¬ 
tering-center  size  0-,  obtained  using  the  average  value  for  the  two  runs, 

2  -8  -3 

£a(0 )/a  «  2.57  x  10  pm  ,  corrected  for  the  longer  wavelength  used  in  the 
IOSA.  The  requirement  for  meeting  the  40  dB  dynamic  range  specification  is 
that  &  be  less  than  about  50  pm.  Our  only  evidence  for  the  value  of  a  is  from 
data  like  that  In  Fig.  19,  which  is  rendered  uncertain  by  the  random  fluctua¬ 
tions  in  the  scattered  intensity  and  by  the  pseudo-scattering  associated  with 
nonuniform  coupling  efficiency.  However,  all  data  of  this  type  is  consistent 
with  values  of  a.  in  the  range  14  pm  -  50  pm,  with  a  *  25  pm  being  an  average 
over  nine  different  runs  using  seven  different  samples.  The  conclusion  de¬ 
rived  from  Fi"  27  is  that  the  IOSA  dynamic  range  will  be  about  42  dB,  or  2  dB 
greater  than  specifications  require.  Although  this  result  is  encouraging,  we 
note  that  the  +  5  dB  error  bars  shown  in  the  figure,  based  on  the  uncertainties 

o 

In  determining  Za(0)/<x  are  critical.  In  the  worst  case  situation,  the  40  dB 
dynamic  range  is  not  met  for  a>12  pm. 

Battelle  Sample  #226 

In  order  to  provide  a  comparison  between  waveguides  fabricated  at 
different  facilities  Battelle  sample  #226  was  subjected  to  the  same  experi¬ 
mental  procedures  that  lead  to  the  results  in  Figs.  23  and  24.  This  waveguide 
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Figure  27.  IOSA  dynamic  range  versus  scattering-center  size 
d,  based  on  the  results  of  Table  IV  for  sample 
Wl. 
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was  formed  using  '''ITS  A  of  Ti  diffused  for  2.5  h  at  965°C  in  moist  flowing 

t  -1  c\ 

oxygen This  fabrication  procedure  is  now  used  routinely  for  most  Battelle 
device-oriented  integrated  optics  programs.  The  surface  quality  of  samples 
fabricated  in  this  way  is  usually  very  good  in  that  the  often  encountered 
orange  peel  texture  of  the  diffused  surface  is  absent.  Sample  226  was  found 
to  be  particularly  good  in  this  regard. 

The  experimental  data  for  relative  scattered  power  versus  scattering 
angle  is  plotted  in  Fig.  28.  For  this  data  an  aperture  was  inserted  in  the 
beam  to  block  the  background  scattered  light.  Accordingly,  the  results 
should  be  valid  over  the  full  angular  range  examined.  We  find  a  clear  in¬ 
crease  in  scattering  associated  with  a  longer  waveguide  path  length  for 
scattering  angles  in  the  range  0.5-2. 2°.  However,  outside  this  interval  the 
data  are  fairly  well  intermixed. 

We  have  analyzed  the  results  to  determine  the  best  Lorentzian  fits  to 

the  data  over  a  specified  angular  range.  Two  different  angular  ranges  were 

selected.  For  0.5° <$<2. 45°  the  data  shows  a  definite  excess  of  scattered 

power  associated  with  a  longer  waveguide  path  length.  For  results  in  this 

2 

interval  the  best  fitting  slopes  of  the  P(0)/AP  versus  <j>  curves  were  used  to 
2  —6  —3 

obtain  Ea(0)/d  *  2.04  x  10  pm  This  is  nearly  two  orders-of-magnitude 

larger  than  the  value  obtained  for  W-l,  indicating  that  IOSA  dynamic  range 

with  the  Battelle  sample  would  be  nearly  20  dB  worse  than  the  value  obtained 

from  Fig.  27  for  the  Westinghouse  sample. 

However,  we  know  of  no  experimental  reason  for  limiting  analysis  to 

angles  less  than  2.45°.  When  we  analyze  all  data  in  the  range  0.5°<<t><6° , 

2 

the  best  fitting  slopes  of  the  (P/Ap)-versus-4>  curves  indicate  a  value 
9  -8  •  q 

£<j(0)/a  ■  4.10  x  10  ym  .  This  is  only  a  factor  1.6  larger  than  we  ob¬ 
tained  with  W-l.  As  with  W-l,  a  value  of  zero  scattering  may  be  consistent 

2 

with  the  data  if  the  true  slopes  of  the  P/AP  versus  <p  curves  are  within 
one  standard  deviation  of  the  statistically  determined  best  values.  It  is 
simply  not  possible  to  state  with  confidence  whether  the  waveguides  are 
good  enough  for  the  IOSA  application,  or  even  whether  the  Westinghouse  fab¬ 
rication  methods  are  to  be  preferred  over  the  Battelle  methods. 
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was  formed  using  ''-ITS  A  of  Ti  diffused  for  2.5  h  at  965°C  in  moist  flowing 
oxygen^  This  fabrication  procedure  is  now  used  routinely  for  most  Battelle 
device-oriented  integrated  optics  programs.  The  surface  quality  of  samples 
fabricated  in  this  way  is  usually  very  good  in  that  the  often  encountered 
orange  peel  texture  of  the  diffused  surface  is  absent.  Sample  226  was  found 
to  be  particularly  good  in  this  regard. 

The  experimental  data  for  relative  scattered  power  versus  scattering 
angle  is  plotted  in  Fig.  28.  For  this  data  an  aperture  was  inserted  in  the 
beam  to  block  the  background  scattered  light.  Accordingly,  the  results 
should  be  valid  over  the  full  angular  range  examined.  We  find  a  clear  in¬ 
crease  in  scattering  associated  with  a  longer  waveguide  path  length  for 
scattering  angles  in  the  range  0.5-2. 2°.  However,  outside  this  interval  the 
data  are  fairly  well  intermixed. 

We  have  analyzed  the  results  to  determine  the  best  Lorentzian  fits  to 

the  data  over  a  specified  angular  range.  Two  different  angular  ranges  were 

selected.  For  0.5° <$<2. 45°  the  data  shows  a  definite  excess  of  scattered 

power  associated  with  a  longer  waveguide  path  length.  For  results  in  this 

2 

interval  the  best  fitting  slopes  of  the  P (0)/AP  versus  <J>  curves  were  used  to 
2  —6  -3 

obtain  Ea(0)/a  =  2.04  x  10  pm  .  This  is  nearly  two  orders-of-magnitude 

larger  than  the  value  obtained  for  W-l,  indicating  that  IOSA  dynamic  range 

with  the  Battelle  sample  would  be  nearly  20  dB  worse  than  the  value  obtained 

from  Fig.  27  for  the  Westinghouse  sample. 

However,  we  know  of  no  experimental  reason  for  limiting  analysis  to 

angles  less  than  2.45°.  When  we  analyze  all  data  in  the  range  O.50<4><6°, 

2 

the  best  fitting  slopes  of  the  (P/AP) -versus-^  curves  indicate  a  value 
2  —8  -3 

lo(0)/<i  =  4.10  x  10  pm  This  is  only  a  factor  1.6  larger  than  we  ob¬ 

tained  with  W-l.  As  with  W-l,  a  value  of  zero  scattering  may  be  consistent 

2 

with  the  data  if  the  true  slopes  of  the  P/AP  versus  curves  are  within 
one  standard  deviation  of  the  statistically  determined  best  values.  It  is 
simply  not  possible  to  state  with  confidence  whether  the  waveguides  are 
good  enough  for  the  IOSA  application,  or  even  whether  the  Westinghouse  fab¬ 
rication  methods  are  to  be  preferred  over  the  Battelle  methods. 


Tigure  28,  Normalized  scattered  power  versus  scattering  angle  for  the  TM0  mode  of  sample  226. 

Dots  correspond  to  14  mm  prism  separation.  Pluses  correspond  to  1  mm  prism  separation 


To  further  illustrate  the  problem,  Fig.  29  shows  data  obtained  using 
Battelle  sample  166,  from  which  the  results  of  Fig.  19  were  obtained  using  a 
different  experimental  technique.  Three  different  prism  separations  were  used 
in  obtaining  the  data  in  Fig.  29,  1  mm,  16  mm,  and  11  mm.  The  last  separation 
was  motivated  by  the  high  scattering  observed  at  16-mm  separation.  It  was 
suspected  that  prism  placement  was  contributing  an  anomalously  large  amount  of 
scattering  to  the  observation.  Judging  from  the  results  of  Fig.  29,  this  was 
apparently  the  case.  However,  the  scattering  at  11-mm  separation  is  anomalously 
low  compared  not  only  to  the  16-mm  result  but  also  to  the  1-mm  result.  For 
this  reason,  we  made  no  attempt  to  statistically  analyze  the  data. 

Surface  Degradation  From  Low-Temperature  Annealing 

In  the  course  of  working  with  the  various  samples  of  the  program,  a 
previously  unobserved  source  of  scattering  came  to  our  attention.  We  had 
noted  early  in  the  program  that  of  the  two  nominally  identical  samples  sent 
to  us  by  Westinghouse,  one  was  of  decidedly  superior  quality  as  judged  by 
in-plane  scattering  levels.  This  is  the  sample  called  W-l  that  we  have  used 
to  take  the  majority  of  data  reported  here.  Late  in  the  program  we  subjected 
the  sample  to  a  series  of  thermal  anneals  to  erase  the  effects  of  optical 
damage  which  may  have  been  generated  in  the  course  of  previous  experiments. 

The  first  anneal  consisted  of  a  15  minute  heat  treat  in  flowing  at  200°C. 

The  sample  was  allowed  to  cool  at  the  natural  rate  of  the  oven,  and  in-plane 
scattering  data  taken  with  1  mm  prism  separation  showed  no  ill  effects  produced 
by  the  anneal. 

Later  a  second  anneal  was  performed  under  identical  conditions  except 
that  the  sample  was  eased  from  the  hot  zone  at  the  rate  of  about  one  inch  per 
ten  minutes.  Figure  30  shows  the  subsequent  intensity  profile  of  the  focused 
beam  after  the  anneal  in  comparison  to  its  previous  appearance. 

After  a  third  anneal  (200°C,  25  min,  flowing  O2)  a  pinging  and  crack¬ 
ling  sound  was  heard  during  withdrawal  conducted  at  the  rate  1  inch  per  15  min. 

In  a  fourth  anneal,  the  withdrawal  rate  was  increased  near  the  end  (30°C  down 
to  room  temperature)  and  the  crackling  sound  was  observed  to  occur  in  coinci¬ 
dence  with  this. 


No  damage  to  the  crystal  surface  was  apparent  to  the  unaided  eye  at 
this  time.  Subsequently  a  fifth  anneal  was  performed  at  500°C  for  1  hour  in 
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Scattering  Angle  in  Air,  Degrees 

Figure  29.  Normalized  scattered  power  versus  scattering  angle  for  the  TM0  mode  of  sample 
166.  +  =  1  mm  prism  separation.  •  =  16  mm  separation.  0  =  11  mm  separation. 


(a) 


(b) 

Figure  30.  Intensity  profile  (linear  scale)  of  focused  beam 
passed  through  sample  Wl.  (a)  Before  anneal 
treatments  (b)  after  two  anneals. 
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flowing  Oj,  followed  by  an  overnight  cool  at  the  natural  rate  of  the  oven.  The 
idea  here  was  that  a  reversible  phase  transformation  could  have  occurred  at 
the  low  temperature  anneals,  and  that  a  longer,  higher  temperature  anneal  might 
restore  the  crystal  to  near  its  original  state. 

Figure  31  shows  the  in-plane  scattered  energy  distribution  for  the 
sample  after  this  heat  treatment  in  comparison  to  previous  results  that  were 
plotted  in  Fig.  24.  It  is  apparent  that  significant  degradation  of  the  wave¬ 
guide  has  occurred  from  the  five  heat  treatments.  Figure  32  shows  a  micro¬ 
graph  of  the  waveguide  surface  taken  at  this  time  that  verifies  this.  Also 
shown  is  a  micrograph  of  the  surface  of  sample  W-2.  It  is  likely  that  same 
phenomena  at  work  in  the  low  temperature  anneals  of  W-l  was  operative  during 
the  fabrication  of  sample  W-2.  The  Implication  is  that  care  must  be  taken  in 
using  the  Westinghouse  fabrication  technique.  Since  Battelle  samples  have  not 
shown  this  problem  and  since  the  main  difference  between  the  Battelle  and 
Westinghouse  fabrication  methods  is  the  use, at  Westinghouse,  of  a  LiNbO^  - 
powder  environment  to  prevent  outdif fusion,  we  speculate  that  this  environ¬ 
ment  is  the  source  of  the  problem. 

Although  we  did  not  subject  Battelle  samples  to  the  same  anneal 
schedule  that  we  used  for  W-l,  a  general  observation  is  that  outdif fusion  pre¬ 
vention  using  flowing  1^0  vapor  leaves  the  waveguide  surface  relatively  free 
of  the  often-observed  orange  peel  texture.  If,  as  we  speculate,  the  presence 
of  orange  peel  correlates  with  surface  strains  that  are  released  in  the  low 
temperature  anneals,  fabrication  methods  which  avoid  orange  peel  are  to  be 
favored.  Unfortunately  we  did  not  perform  a  careful  micrographic  examination 
of  samples  W-l  and  W-2  when  they  were  received  to  see  if  they  had  orange  peel. 
Because  of  this  and  because  we  did  not  subject  sample  226  to  the  same  anneal 
schedule  as  W-l,  we  do  not  have  a  sufficient  set  of  evidence  from  which  to  draw 
definite  conclusions,  but  we  believe  that  our  speculations  warrant  attention. 

For  completeness  we  note  that  the  data  in  Fig.  28  for  Battelle  sample  226 
was  obtained  following  a  low- temperature  anneal,  but  it  is  not  known  whether 
the  sample  was  cooled  rapidly  enough  to  provide  an  adequate  test  for  thermal 
stability. 
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Figure  32.  Micrographs  of  the  surface  of  the  Westinghouse  samples 

(a)  Sample  W-l,  after  five  low- temperature  anneals. 

(b)  Sample  W-2,  after  one  low-temperature  anneal. 
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OTHER  IN-PLANE  SCATTERING  EXPERIMENTS 


The  Use  of  End-Fire  Input  and  Edge  Output  Coupling 

Other  experiments  to  measure  in-plane  waveguide  scattering  were  per¬ 
formed  using  different  geometries  to  circumvent  the  problem  of  scattering  from 
prism  coupling.  The  obvious  modification  of  this  type  involves  the  use  of  end-fire 
coupling  to  replace  prism  coupling.  However,  imperfect  polishing  of  the  edges 
leads  to  scattering  that  is  in  its  effect  no  different  than  that  associated 
with  prism  coupling.  If  the  edge  has  an  rms  roughness  a  and  an  autocorrela¬ 
tion  length  a,  an  analysis  similar  to  that  used  for  the  derivation  of  Eq.  38 
shows  that  the  associated  scattered-power  distribution  is  given  by 

AP/P(0)  -  (16/tt)  k2(n-l)2a2(a/W)L(k  a*).  (64) 

o  o 

For  k  *  (2ir/0.83)  pm  \  a  =  20  pm,  W  =  4  mm,  and  n  =  2.2,  an  on-axis  scatter- 

°  O 

ing  level  of  30  dB  is  consistent  with  a  surface  roughness  of  a  *  2.62  A.  This 
is  much  rougher  than  the  levels  that  may  be  achieved  in  polishing  flat  surfaces 
of  LiNbO^1^ However,  the  special  problems  encountered  in  polishing  sharp  edges 
on  the  material  may  lead  to  roughness  values  of  this  magnitude. 

We  have  performed  the  experimental  comparison  in  which  a  focused 
beam  was  end-fired  into  Westinghouse  sample  W1  and  output  coupled  using  either 
a  prism  or  diffraction  from  the  polished  waveguide  termination  edge.  The 
experimental  geometry  is  shown  in  Fig.  33  for  the  case  of  prism  output  coupl¬ 
ing.  The  results  of  the  comparison  experiment  are  displayed  in  Fig.  34.  It 
is  evident  that  scattering  from  polished-edge  output  coupling  is  greater  than 
that  from  prism  output  coupling.  As  a  result  of  this  measurement,  we  did  not 
further  investigate  the  former  output-coupling  method.  The  theoretical  curve 
for  the  case  of  prism  output  coupling  will  be  discussed  in  the  following 
paragraphs . 
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Distance  from  Focus,  cm 

Effect  of  output  coupling  method  on  scattered  energy 
distribution.  #,  End-fire  input  and  edge  output  cou 
pling.  o,  End-fire  input  and  prism  output  coupling. 
Experimental  geometry  is  that  of  Figure  33. 


The  Image  Geometry  for  Studies  of  In-Plane  Scattering 

From  Fig.  33,  it  may  be  noted  that  the  focused  beam  incident  on  the 
detector  is  really  an  image  of  the  focused  beam  coupled  into  the  waveguide. 

This  is  distinguished  from  the  conventional  measurement  using  a  collimated 
input  beam,  in  which  case  the  focused  beam  at  the  detector  is  a  Fourier  trans¬ 
form  of  the  input  beam. 

The  significance  of  this  is  illustrated  in  Fig.  35.  in  the  Fourier- 
transform  geometry,  all  rays  scattered  at  the  same  angle  are  projected  to  the 
same  off-axis  position  y  in  the  focal  plane.  Because  of  this  feature,  the 
shape  of  the  scattering  cross  section  is  known  directly  from  the  shape  of 
the  scattered-energy  distribution.  This  feature  is  absent  in  the  imaging 
geometry,  because  the  position  y  to  which  a  given  scattered  ray  is  projected 
depends  not  only  the  angle  of  scattering  but  also  on  the  location  of  the  scat¬ 
tering  center  in  the  waveguide.  This  makes  it  necessary  to  do  a  bit  of  model¬ 
ing  to  extract  scattering  parameters  from  the  data;  however,  an  advantage  is 
that  rays  scattered  from  the  point  of  input  coupling  are  imaged  to  y  *»  0  and 
do  not  contribute  to  the  scattered  m-line  distribution.  This  enables  us  to 
eliminate  background  scattering  from  input  coupling  from  our  measurement,  and 
is  very  useful  if  output  coupling  can  be  accomplished  with  little  additional 
scattering  added  to  that  from  the  waveguide.  Assuming  that  this  is  a  possi¬ 
bility,  we  made  a  careful  analysis  of  those  results  in  Fig.  34  obtained  using 
prism  output  coupling,  and  we  also  employed  the  image  geometry  to  study  other 
waveguide  samples  as  well. 

Interpretation  of  Image-Geometry  Results 

The  analysis  to  determine  scattering  center  parameters  from  image- 
geometry  results  is  similar  to  that  employed  to  calculate  the  scattered  signal 
from  che  regions  of  the  IOSA  where  the  guided  beam  is  either  converging  or 
diverging  (regions  A  and  D  of  Fig.  1).  The  approach  is  to  divide  the  beam 
path  in  the  waveguide  into  increments  of  area  and  to  calculate  the  scattered 
power  from  each  increment  received  at  a  selected  pixel  in  the  detector  array. 
The  total  power  is  calculated  by  integrating  over  all  increments  of  area  in 
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Figure  35.  (a)  Conventional  Fourier  -  transform  approach  for  measuring 

in-plane  scattering,  in  which  all  rays  scattered  at  angle  9 
are  projected  to  the  same  point  on  the  m-line.  (b)  Imaging 
method  for  measuring  in-plane  scattering,  in  which  rays  scat¬ 
tered  at  the  point  of  input  coupling  (Z  *  0)  do  not  contrib¬ 
ute  to  the  measured  m-line. 


the  beam  path.  If  the  pixel  has  an  aperture  Ay  and  is  a  postive  distance  y 
from  the  optical  axis,  the  received  power  in  the  case  of  RGD  scattering  centers 
is 

AP  -  PoIa(0)M_1  [y(I+-I_)  +  (Ay/2)  (I+  +  I_) ] 

I+  =  iDjta/LXy  +  Ay/2)  /M] 

I(u)  *  u  ^  tan  1u  +  £n(l  +  u  2)^2.  (65) 

In  this  expression  Pq  is  the  total  power  and  M  *  Si/SQ  *s  c^e  magnification 
associated  with  the  imaging  geometry,  approximately  equal  to  25  in  our 
experiments . 

For  all  values  of  y  that  we  measured,  an  excellent  approximation 
to  Eq .  65  is 


AP  =  PQ  Io(0)  M  1Ay  «n[l  +  (ML/kfly)2]1/2  (66) 

The  ratio  of  scattered  power  to  the  peak,  on-axis  unscattered  power  is 

AP/P(0)  *  (2/ir) 1/2  (FAq/W)  Io(0)  £n[l  +  (ML/kfly) 2]1/2  (67) 

The  theoretical  curve  in  Fig.  34  is  based  on  the  parameter  values 

Eo(0)  ■  0.51  mm  *  and  a  *  16.5  pm.  The  range  of  data  acquired  is  sufficiently 

large  that  values  for  both  parameters  can  be  obtained  separately, 

2 

rather  than  in  the  combination  Eo(0)/a  ,  as  was  the  case  with  the  Fourier- 
transform  geometry.  However,  when  we  calculate  the  ratio  using  the  above 
values  for  la(0)  and  a,  we  obtain  1.87  x  10  ^  ym  ^ .  This  is  about  two  orders 
of  magnitude  greater  than  we  achieved  using  the  Fourier  transform  geometry. 

We  were  never  able  to  resolve  why  the  scattering  should  be  so  much  greater  for 
the  image  geometry,  and  we  abandoned  it  as  a  measurement  technique  later  in 
the  program. 
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IOSA  DYNAMIC-RANGE  FORECAST 


The  IOSA  dynamic  range  specifications  are  prescribed  at  the  fourth 
off-axis  pixel  on  either  side  of  a  central  illuminated  pixel.  This  has  been 
taken  into  account  in  generating  previous  figures  in  which  dynamic  range  has 
been  plotted  against  various  scattering  parameters.  It  is  important,  however 
to  know  how  the  dynamic  range  behaves  as  a  function  of  pixel  number,  given 
the  scattering  parameters  measured  in  this  program.  The  information  is  pro¬ 
vided  here.  We  emphasize  that  the  calculations  are  only  as  good  as  the 
measurements,  and  the  measurements  were  plagued  with  experimental  uncertainties 
from  background  scattering  that  were  not  successfully  resolved.  Consequently 
the  results  are  best  estimates  that  may  or  may  not  hold  up  as  spectrum  ana¬ 
lyzers  are  actually  made  and  tested. 

The  relation  to  be  plotted  is  10  log  [AP(m)/P]  versus  m,  where 
APtQt(m)  is  the  total  scattered  power  received  at  pixel  m  and  P  is  the  power 
launched  into  the  waveguide  mode.  The  ratio  AP(m)/P  is  calculated  from  IOSA 
design  and  scattering  parameters  according  to  Eq.  25.  Values  for  these  are 

listed  in  Table  V.  The  design  parameters  are  relevant  to  the  Westinghouse 
(4) 

IOSA  effort,  while  the  scattering  parameters  are  those  obtained  for 
sample  W-l. 

With  regard  to  specification  of  the  scattering  parameters,  a  word 

2 

of  explanation  is  in  order.  For  two  different  runs  we  found  la(0)/n  ■ 

(6.2  +  2.4)  x  10  ®  ym  ^  and  (-1.0  +  1.9)  x  10  ®  ym-^.  Assigning  equal  weight 

—8  —3 

to  each  run,  the  average  is  2.6  x  10  ym  .  The  standard  deviation  calcu- 

—8  —3 

lated  from  the  input  standard  deviations  is  2.2  x  10  ym  while  that  calcu- 

_g  _  q 

lated  from  the  two  point  sample  population  is  5.0  x  10  ym  .  If  we  base  our 

2  -8  -3 

calculations  of  IOSA  dynamic  range  on  the  average  Eo(0)/d  *  2.6  x  10  ym  , 

there  is  a  reasonable  chance  that  our  calculated  dynamic  range  values  are  too 
high  by  3  to  5  dB,  depending  on  which  value  for  standard  deviation  we  select. 
Equally  consistent  with  our  results  is  the  possibility  that  our  calculated 
values  are  too low  by  the  same  amount.  The  possibility  that  the  waveguide 
scattering  is  acutally  zero  is  also  consistent  with  our  results;  however, 
the  preponderance  of  our  data  taken  over  three  years  shows  an  average  scatter¬ 
ing  greater  than  zero.  Therefore,  we  are  inclined  to  minimize  this  possi¬ 
bility.  Error  bars  of  +  5  dB  shown  in  the  figures  are  based  on  these  delibera¬ 
tions. 


TABLE  V.  DESIGN  AND  SCATTERING  PARAMETERS  USED  TO 
CALCULATE  SCATTERING  LIMITATIONS  TO  IOSA 
DYNAMIC-RANGE 


Lens  Focal  Length 

F 

27 . 2  mm 

Lens-Lens  Separation 

L 

18.0  mm 

Detector  Aperture  and  Spacing 

<SY 

12.0  pm 

Diffraction  Angle 

* 

0.066  radians 

Diffraction  Efficiency 

n 

0.05 

Free-Space  Wavelength 

Xo 

0 . 83  m 

Refractive  Index 

n 

2.2 

Scattering  Parameter 

Eo(0)/a2 

1.15  x  10~8  pm-2 
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In  Table  V  the  value  for  Za(0)/a  have  been  corrected  for  the  fact 
that  the  IOSA  will  operate  with  light  of  0.83  pm  wavelength.  Our  experiments 
were  performed  at  0.633  po,  and  for  RGD  scattering,  o(0)  is  proportional  to 

^  as  described  by  Eq.  39.  Therefore  Zo(0)/a^  ■  2.6  x  10  ®  x  (0.633/0.83)^  * 
—8  —3 

1.15  x  10  pm  .  This  value  is  used  in  Table  V.  Note  that  in  general  there 
will  be  a  wavelength  dependence  of  the  index  parameter  Sn  which  is  an  impor- 

O 

tant  factor  in  Eq.  39.  However,  we  have  neglected  this  possibility  in  our 
analysis.  The  results  of  our  analysis  are  shown  in  Fig.  36,  where  dynamic 
range  versus  pixel  number  is  plotted  with  scattering  center  size  as  a  parameter. 
The  range  of  scattering  center  size  selected,  from  10  to  50  pm,  is  consistent 
with  our  experimental  evidence.  Interestingly,  dynamic  range  is  not  sensitive 
to  this  parameter  at  large  pixel  numbers.  Also  note  the  slight  asymmetry  in 
the  curves  associated  with  the  contribution  to  scattering  from  the  0-order 
undiffracted  beam.  This  beam  is  focused  at  a  position  corresponding  to  pixel 
position  150.  Owing  to  the  rapid  fall  off  of  scattered  intensity  with  scat¬ 
tering  angle  for  RGD  scattering  centers,  the  contribution  from  this  source  is 
not  significant  even  though  it  carries  95%  of  the  total  unscattered  IOSA  beam 
power. 

These  results  suggest  that  the  IOSA  dynamic  range  specifications  are 
met  by  sample  W-l  and  the  Westinghouse  design.  Even  in  the  worst  case  repre¬ 
sented  by  the  +  5  dB  error  bar  in  Figure  36,  40  dB  dynamic  range  is  achieved 
at  the  eighth  off-axis  pixel. 
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Figure  36.  IOSA  dynamic  range  versus  pixel  number  with  scattering  center  size  as  a  parameter. 

The  curves  are  based  on  the  Westinghouse  design  and  waveguide  scattering  results 
shown  in  Table  V.  An  error  bar  indicates  the  amount  each  curve  may  be  raised  or 
lowered  to  account  for  experimental  uncertainties  in  the  measurement  of  E o(0)/a^. 


V.  SUMMARY  AND  CONCLUSIONS 


This  report  has  described  the  results  of  a  research  program  to 
characterize  and  reduce  the  effects  of  scattering  in  Ti-diffused  LiNbO^  wave¬ 
guides.  This  research  has  been  performed  for  the  Naval  Research  Laboratories 
to  extend  the  results  of  two  previous  programs  on  the  same  subject  conducted 
for  the  Air  Force.  The  extent  of  government  support  shows  the  perceived 

importance  of  scattering  to  the  success  of  planar  Integrated  optics  as  a 
signal-processing  technology.  All  who  have  been  involved  with  the  technology 
and  the  problem  of  scattering  have  appreciated  its  difficulty.  We  neverthe¬ 
less  expected  that  the  results  of  this  work,  as  exemplified  by  Fig.  36,  would 
be  more  definitive  if  not  more  optimistic.  Because  of  experimental  problems, 
the  program  emphasized  the  measurement  of  waveguide  scattering  rather  than 
the  understanding  of  its  causes  and  the  means  for  its  reduction.  In  the  future, 
progress  in  these  areas  will  come  from  those  who,  viewing  scattering  from  a 
fresh  perspective,  discover  new  experimental  approaches  to  the  problem.  We 
take  considerable  satisfaction  in  providing  these  workers  with  a  sound  founda¬ 
tion  on  which  to  base  their  research. 

First  of  all,  our  evidence  suggests  the  dominance  of  Rayleigh-Gans- 

Debye  scattering  associated  with  scattering  centers  in  the  size  range  10-50  pm 

and  typical  effective  refractive  index  inhomogeneity  10  .  These  numbers 

alone  suggest  the  difficulty  in  directly  observing  the  scattering  centers 

using  diagnostic  tools.  Our  best  assessment  is  that  the  scattering  sources 

lie  close  to  the  waveguide  surface  and  that  the  low  value  of  6n^  results  from 

their  small  depth  normal  to  the  waveguide  plane.  Two  candidate  mechanisms  for 

such  scattering  centers  are  surface  roughness  and  Li-Ti-0  compound  formation 

during  diffusion.  Stylus  profilometers  are  adequate  for  the  measurement  of 

surface  roughness  and  this  is  certainly  one  diagnostic  approach  recommended 

for  future  scattering  research.  SIMS  analysis  has  been  useful  for  verifying 

(2) 

the  occurence  of  anomalously  high  T1  levels  at  the  waveguide  surface,  but 
this  technique  is  not  likely  to  have  the  lateral  resolution  necessary  to 
identify  and  study  individual  scattering  sites.  Nevertheless,  by  combining 
stylus  measurements,  SIMS  measurements  and  in-plane  scattering  measurements 
it  is  believed  that  valid  conclusions  about  the  dominant  source  of  scattering 
can  be  drawn.  This  was  the  plan  of  the  current  program,  and  it  is  still 
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recommended  for  future  work  If  the  experimental  problems  in  generating  good 
in-plane  scattering  data  may  be  solved. 

Assuming  that  surface-localized  scattering  centers  are  dominant  in 

Ti-LINbOj  waveguides,  scattering  reduction  is  simply  achieved  by  reducing  the 

waveguide  modal  field  value  at  the  surface.  This  generally  means  that  longer, 

hotter  diffusion  treatments  produce  better  waveguides,  other  things  being 

equal.  But  other  things  are  generally  not  equal.  Diffusion  time  and  tempera- 

(14) 

ture  have  an  effect  on  surface  compound  formation,  while  surface-compound 
formation  appears  to  have  an  effect  on  surface  roughness. Thus  the 
interplay  of  phenomena  is  complex.  Each  datum  point  to  be  used  in  analyzing 
the  situation  requires  the  fabrication  of  a  new  waveguide  and  the  subsequent 
performance  of  stylus,  SIMS,  and  in-plane  scattering  measurements,  for  example. 
Thus,  the  approach  to  the  problem  is  likely  to  be  expensive  and  time  consuming. 
Based  on  the  limited  effort  of  this  program,  we  have  not  noticed  significant 
differences  between  waveguides  to  result  from  variation  of  fabrication  pro¬ 
cedures,  within  reasonable  and  commonly  accepted  standards.  One  Westinghouse 
waveguide  was  superior  to  Battelle  waveguides,  but  degraded  as  the  result  of 
a  low  temperature  thermal  anneal.  Another  Westinghouse  waveguide  was  inferior 
to  Battelle  samples.  In  short,  there  has  been  little  to  recommend  one  fabrica¬ 
tion  technique  over  another. 

Despite  this  unen thus las  tic  summary  of  our  current  understanding  of 
the  scattering  problem,  one  encouraging  result  stands  out:  when  our  best 
measurements  are  made  and  analyzed  to  the  best  of  our  ability,  the  scattering 
limitations  to  IOSA  dynamic  range  are  forecast  to  be  within  specifications.  By 
doing  nothing  more  to  understand  and  reduce  waveguide  scattering,  we  may  still 
produce  waveguides  good  enough  for  many  integrated-optics  applications.  And 
they  will  undoubtedly  get  better  as  people  make  and  use  more  waveguides  for 
an 'expanding  set  of  applications. 
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